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Abstract 

In  separate  parts,  this  report  addresses:  (I)  the  use  of  multiscale  modulation  to  obtain  bandwidth-, 
efficient  modulation,  and  the  performance  of  such  modulation  in  fading  channels;  (2)  The  analysis  of 
transmission  through  fading  channels,  and  approaches  to  the  design  of  transmitter  and  receiver  waveforms 
which  will  perform  optimally  in  fading  channels.  In  the  first  case,  it  is  shown  that  using  wavelet  and/or 
scaling  waveforms  can  lead  to  modulation  which  is  comparable  to,  or  better  than,  many  conventional 
schemes,  including  GMSK.  Performance  in  fading  channels  is  similar  to  that  of  other  methods.  The 
performance  of  raultiscale  signaling  is  also  addressed  using  the  matched-filter  bound  method.  In  the 
second  part,  the  signal  design  problem  is  formulated  as  an  optimization  problem  and  an  analysis  approach 
outlined  which  leads  to  maximization  of  SNR  in  fading  channels.  Inclusion  of  constraints  for  good  spectral 
localization  is  also  possible. 


1  Introduction 

In  this  report,  results  from  two  related  directions  of  research  axe  presented.  First,  results  related  to  wavelet 
modulation  are  discussed.  Second,  results  related  to  waveform  design  for  communication  over  fading  channels 
are  presented. 

Wavelet  modulation  employs  as  its  baseband  basis  waveforms  the  shift-  and  scale-orthogonal  scaling 
functions  and  wavelets,  typified  by  Daubechies  wavelets  (see,  e.g.,  [1]).  Because  of  the  orthogonal  properties, 
it  is  possible  to  employ  multiscale  modulation,  sending  information  on  multiple  time  scales. 

One  of  the  potential  benefits  of  wavelet  modulation  is  that,  as  basis  functions  with  longer  support  are 
used  (  resulting  in  multiple  time-overlapped  waveforms),  the  smoothness  of  these  waveforms  leads  to  better 
spec t nil  localization,  as  measured  by  the  fraction  of  out-of-band  power  (FOOBP),  than  many  other  popular 
signal  formats,  including  GMSK. 

As  one  of  the  research  hypothesis,  it  was  conjectured  that  transmitting  on  multiple  scales  and  over  long 
time  frajn<*s  might  yield  a  performance  benefit  in  channels  experiencing  fast  fading.  The  conjecture  arose 
with  the  idea  that  a  short  fading  interval  might  affect  only  a  small  fraction  of  the  overall  waveform.  A 
significant  part  of  the  research  project  was  devoted  to  exploring  this  question.  The  findings  are  that  in 
some  cases,  there  may  be  a  modest  benefit  in  a  fading  channel,  in  comparison  to  conventional  (flat-pulse) 
signaling  However,  the  improvements  seemed  to  apply  to  such  a  narrow  set  of  specific  cases,  and  in  any 
event  he  sufficiently  modest,  that  it  seemed  difficult  justify  the  use  of  these  basis  functions  on  that  criterion 
alone. 

This  leads  to  the  question  of  whether  there  might  be  other  waveforms  which  could  be  specifically  designed 
to  be  effective  in  fading  channels  and  which  might  be,  furthermore,  constrained  to  have  a  sufficiently  high 
degree  of  spectral  localization  that  they  could  be  used  in  practical  channels.  This  therefore,  forms  the  turning 
point  for  the  second  part  of  the  research,  into  waveform  design  for  fading  channels.  As  of  the  time  of  this 
report,  approaches  to  the  problem  have  been  formulated,  as  will  be  discussed,  but  the  performance  of  these 
methods  compared  with  other  channels  remains  to  be  accomplished. 

The  organization  of  this  report  is  as  follows.  Sections  2  through  3.2  present  the  basic  notation,  spectral 
efficiency,  and  analysis  for  communication  over  fading  channels  using  differential  detection.  In  section  4, 
the  matched  filter  bound  is  employed  to  analyze  multiscale  transmission  over  several  fading  profiles.  It 
is  found  that  the  performance  is  somewhat  better  than  raised  cosine  waveforms  in  many  cases,  but  again 
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the  improvement  is  at  best  moderate.  Portions  of  these  materials  (greatly  abbreviated)  were  presented  at 
the  International  Telemetering  Conference  (ITC)  [2].  In  section  6,  the  concept  of  adapting  the  transmitted 
spectrum  by  tuning  wavelet  parameters  is  briefly  introduced. 

In  Part  II,  we  consider  several  aspects  of  the  analysis  and  design  of  signal  waveforms  for  fading  channels. 
Section  7  provides  a  discussion  of  modeling  for  transmission  over  fading  channels,  describes  some  fading 
channel  models,  and  provides  a  literature  survey  or  related  work.  This  work  has  not  yet  appeared  in  print. 
In  section  8,  an  analysis  of  transmission  through  a  multiplicative  fading  channel  is  performed  using  two 
methods  of  analysis.  This  section  is  strictly  a  matter  of  analysis,  undertaken  to  understand  more  clearly  the 
nature  of  fading.  This  section  has  not  been  published  elsewhere.  In  section  9,  a  method  is  introduced  to 
design  signal  waveforms  which  are  optimal  with  respect  to  fading  without  the  assumption  of  synchronization. 
This  design  is  extended  in  section  10  to  the  question  of  designing  the  signaling  waveform  and  the  receiver 
waveform  to  achieve  maximum  signal  to  noise  ratio  when  perfect  synchronization  is  assumed.  The  design 
methodology  is  outlined  but,  due  to  the  fact  that  research  is  still  underway,  waveforms  and  performance 
results  are  not  presented. 

A  set  of  appendixes  describes  some  details  of  computations  need  in  sections  3.1  and  4. 

In  summary,  this  research  summarized  in  this  report  offers  the  following  contributions: 

1 .  A  model  for  using  wavelet  waveforms  and  employing  their  shift  and  scale  orthogonality  properties  to 
obtain  flexible,  bandwidth-efficient  modulation  waveforms. 

2.  An  analysis  of  the  performance  of  these  signals  in  Rician  fast  fading  channels  with  differential  detection, 
showing  that  the  waveforms  perform  well,  but  with  improvements  that  are  only  moderate. 

3.  A  matched-filter  bound  analysis  of  multiscale  transmission  through  some  common  channel  models, 
showing  performance  comparable  (and  usually  somewhat  better)  than  a  similar  analysis  for  raised 
cosine  signal  waveforms. 

4.  An  introduction  to  the  concept  that  the  waveforms  can  be  modified  parametrically  to  match  the 
waveform  to  the  channel  spectrum. 

5.  An  analysis  of  a  multi-chip  fading  channel  model. 

G.  A  design  approach  for  finding  the  optimal  transmission  waveform  and  receiver  filter  matched  to  a 
channel  with  a  prescribed  fading  autocorrelation  function. 
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Part  I 

Wavelet  Modulation 

2  Notation 

A  notation  is  established  here  which  is  used  throughout  this  first  part.  A  multiresolution  analysis  [3,  1,  4,  5] 
is  formed  by  a  sequence  of  closed  sub-spaces  Vj  of  L2(K)  which  are  nested  according  to 

•  •  ■  V2  C  Vi  C  V'0  C  V-i  C  V_2  •  •  •  . 

The  spaces  have  the  following  pertinent  properties  (among  others):  the  scale  invariance  property, 

/(f)  G  Vm  /(2mt)  G  Vo 

the  shift  invariance  property, 

/(f)  G  Vo  =>  /(f  -  n)  €  V0  for  all  n  G  Z, 

and  the  existence  of  a  basis  function  0(f)  such  that 

Vo  =  span {<j>(t  -  n);  n  €  Z} 

and  where  span{  }  is  the  linear  span  of  the  set  of  functions  in  the  argument.  The  basis  functions  additionally 
satisfy  the  shift  orthogonality  property 

{4>(t  -  j),<f>(t  -  k))  =6jik, 

where  the  inner  product  is  (/,  g)  =  f  f(t)g(t)  dt  and  Sj,k  is  the  usual  Kronecker  delta  with  integer  arguments. 
The  function  <t>{t)  is  known  as  a  scaling  function  and  can  be  shown  to  be  a  lowpass  function.  The  shift 
orthogonality  property  of  the  scaling  function  (1)  can  be  shown  to  be  equivalent  to  the  Fourier  domain 

constraint 

°o 

y  \(p(u )  +  2Attt)|2  =  1,  (2) 


where  oU')  is  the  Fourier  transform  of  0(f)  [1].  The  constraint  (2)  is  equivalent  to  the  Nyquist  zero  ISI 
criterion  (sc<\  e.g.  [6,  p.  561]),  so  the  class  of  scaling  functions  includes  Nyquist  signals,  including  the 

square* root  raised- cosine  waveform. 

Let  Wj  be  the  orthogonal  complement  of  Vj  in  Vj- 1;  i.e., 

Vj-i  =  VjQWj,  (3) 

where  rv  indicates  the  direct  sum.  It  follows  that  Wj  is  orthogonal  to  Wj >  for  j  /  j‘  and  that  for  j  <  J, 

J-j-i 

Vj  =  Vj  ©  0  Wj-k  (4) 

k-0 

where  all  the  subspaces  on  the  right  are  orthogonal  to  each  other. 

It  can  be  shown  that  there  is  a  function  0(f)  such  that 

Wo  =  span{0(f  -  n);n  €  Z}. 

The  function  v(f )  is  a  wavelet  function  and  is  a  bandpass  signal.  From  the  properties  of  the  multiresolution 
analysis,  it  follows  that  0(f)  G  V0  C  V-i  must  be  a  linear  combination  of  shifts  of  0(2f)  €  V_i, 


;  v  - 1 

0(f)  =  ^2  ck<t>(2t  -  k) 
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(6) 


and  ip(t)  €  Wo  CV-i  can  be  expressed  as  the  linear  combination 

N-l 

jp(t)  =  dk<f>(2t  -  k) 
k^o 

for  appropriately  selected  sequences  {c*}  and  {d*}.  Let  —  2j^2<f>(2H  —  k)  and  ipj,k{t)  =  2j/2ip(2jt  -  k). 

Then  by  the  properties  of  the  multiresolution  analysis, 

(^fcW.lh.rnW)  =  0  3  >  1  and  j, k,l,m  €  Z 

and 

In  (5),  if  the  number  of  coefficients  N  is  finite  then  the  function  has  compact  support.  While  compact 
support  is  not  necessary  for  communication  applications  (e.g.  the  square-root  raised-cosine  does  not  have 
compact  support),  compact  support  does  simplify  some  aspects  of  implementation  and  analysis.  It  is  assumed 
throughout  this  report  that  the  scaling  function  employed  has  support  on  f  E  [0,g].  Compactly-supported 
scaling  functions  give  rise  to  compactly-supported  wavelet  functions.  Wavelet  and  scaling  functions  param¬ 
eterized  by  N  coefficients  in  (5)  and  (6)  are  said  to  be  in  the  Djv  family,  where  D  represents  Daubechies  [1] 
and  N  must  be  even.  The  support  and  smoothness  of  the  <j>  and  ip  functions  increase  with  N.  Observe  that 
the  wavelet  function  with  two  coefficients  is  the  familiar  Haar  function. 

Let 

<t>' t (t)  =  —/=<t>(t/T) 

and 

t (t)  =  -j=ii>{t/T). 

These  scaled  functions  can  be  used  to  transmit  one  symbol  of  information  every  T  seconds  according  to  the 
baseband  signal  model 

Sd(t)  =  y/2Sg(t),  (7) 

whore 

ffW  =  -  *T),  (8) 

i€Z 

and  S  is  the  signal  energy,  a*  €  C ,  and  C  is  some  real  or  complex  signal  constellation  with  2nh  signal  points. 
In  (8),  T  is  said  to  be  the  symbol  time,  being  the  time  interval  between  transmission  of  adjacent  symbols, 
even  though  the  duration  of  the  symbol  is  qT  seconds.  Each  scaling  function  transmits  a  single  symbol.  By 
the  shift  orthogonality  property,  the  outputs  of  a  matched  filter  provide  sufficient  information  to  detect  the 
received  signal  independently  of  any  other  symbols.  This  is  similar,  in  some  respects,  to  partial  response 
signaling,  in  which  the  signal  is  spread  out  over  several  symbol  times.  However,  unlike  partial  response 
signaling,  the  detector  does  not  need  to  traverse  a  trellis  or  keep  track  of  previous  symbols,  because  of  the 
orthogonality  between  signals.  Like  a  Nyquist  zero  IS1  signed,  the  sampled  outputs  of  the  matched  filter  are 
orthogonal  to  signals  at  other  symbol  times  and  on  other  scales. 

Using  the  multiresolution  property  of  (3),  a  signal  spanning  the  same  space  can  be  written  using  basis 
functions  far  and  V2 r,  catch  with  independent  data  streams: 

g(t)  =  r{t  -  2 it)  +  -  2iT).  (9) 

t€Z 

If  the  signal  amplitudes  and  a ^  are  each  drawn  from  constellations  containing  2Ub  signals,  then  the 
data  rate  for  the  signal  of  (9)  is  the  same  as  for  the  signal  of  (8),  even  though  the  symbol  times  are  twice 


V 


as  long.  The  multiresolution  property  can  be  applied  recursively  Ns  times  to  the  scaling  function  in  (9),  as 
suggested  by  (4).  The  transmission  using  Ns  scales  can  be  written  as 

g(t)  =  ^  —  > 2aT)  +  ^  a<M<&2*.r(f  —  i2N‘T).  (10) 

<r=l  t€Z  »ez 

Observe  that  if  the  transmission  on  each  scale  is  binary,  then  the  energy  per  bit  is  the  same  for  bits  on  all 
scales.  This  means  that  the  power  per  bit  decreases  with  increasing  scale  length. 

For  the  sake  of  the  computations  below,  it  will  be  convenient  to  introduce  more  compact  notation.  Let 


j  ih -tW  1  <  <  K, 

°  +  1 

be  used  to  represent  the  scaling  function  and  the  wavelets  on  the  different  scales,  depending  on  the  value  of 
a.  Throughout  the  report,  let  &’  =  min(cr,  A » )  and  s'  =  min(s,  Nt)  denote  the  actual  scale  number  and  let 
Ta,  =  2 °'T.  The  signal  amplitudes  on  scale  a  at  the  ith  symbol  be  indexed  by  multiples  of  T  ,  with 

,  _  /  1  <  o  <  Ns 

W*  ~  \  a*,t-  a  =  N,  +  l. 

It  is  also  convenient  to  let  qa  =  2 ”q  (the  length  of  the  support  of  the  signal  on  scale  a),  with  qa  =  qa‘  when 
a  =  .V,  +  1.  Using  these  notational  conventions,  the  signal  (10)  can  be  written  as 


JV.  +  l 

5(f)  =  ^  ^<7.2^'  (f 

<r  =  1  i€Z 


iT„.) 


The  signal  in  (11)  can  also  be  written  as 


N.t- 1 

g(t)  =  ^  ^  ~  bff'iPutf  —  iT) 

<7=\  i£2<r'  1 


(11) 


(12) 


Figure  2(a)  illustrates  the  time-overlap  of  multiscale  signals  for  the  particular  case  when  N,  -  2  and  q  =  5, 
using  the  notation  of  ( 12).  At  every  instant  in  time  (AT,  +  1  )q  symbols  are  being  transmitted.  The  figure  also 
illustrates  the  bit  indexing  convention.  Figure  3  shows  another  point  of  view  for  multiscale  transmission. 
Rather  than  focusing  on  the  overlap  across  the  different  signals,  it  illustrates  the  intervals  in  which  new  signals 
are  launched.  The  overall  signal  pattern  in  figure  3  repeats  every  Ns  symbol  times.  The  group  of  signals 
within  a  repeating  pattern  is  called  a  supersymbol ,  with  length  N,T  seconds  and  containing  2 ;V*  individual 
symbols.  The  constituent  symbols  in  a  supersymbol  may  be  considered  as  components  in  a  multidimensional 
signal  space.  By  sending  information  on  a  subset  of  the  full  set  of  scales,  a  variety  of  signal  dimensionalities 
may  be  obtained. 

It  is  interesting  to  examine  the  spectral  efficiency  of  multiscale  signaling.  Due  to  the  multiresolution 
property  of  the  scaling  functions,  it  can  be  shown  [7]  that  if  the  signed  constellation  on  each  scale  is  the 
same,  the  power  spectrum  of  the  multiscale  signal  (11)  is  the  same  as  the  power  spectrum  of  the  single-scale 
signal  (8).  The  spectral  efficiency  of  the  signal  may  be  determined  by  means  of  the  fractional  out-of-band 
power  (FOOBP).  defined  by 


77(6)  =  1  - 


f”'%  Gl/ldf 

IZeOmdf 


where  G(f)  is  the  power  spectral  density  of  the  signed, 


G(f)  = 


m*n\2 

T 


5 


and  <£(2*7)  is  the  Fourier  transform  of  4>{t).  Figure  4  illustrates  the  FOOBP  for  transmission  using  scaling 
functions  from  the  D4,D6,D8,  and  D10  families.  Also  shown,  for  comparison,  is  the  FOOBP  for  flat- 
topped  pulse  signaling,  MSK,  SFSK,  and  GMSK  with  h  =  0.36  [8].  For  each  waveform,  4-PSK  signaling  (2 
bits/symbol)  is  used.  Being  a  continuous  signal,  the  Z)4  scaling  function  is  more  spectrally  localized  than 
the  flat-topped  pulsed,  but  not  as  much  as  MSK.  Scaling  functions  having  more  coefficients,  hence  better 
regularity,  have  better  localization.  The  D6  scaling  function  has  spectral  localization  similar  to  that  of  MSK, 
and  £>s  and  Dl0  do  better  than  MSK.  The  longer  scaling  functions  do  as  well  as  GMSK  or  SFSK  or  other 
spectrally  efficient  techniques.  Even  longer  scaling  functions  are  possible,  resulting  in  even  better  spectral 
localization.  The  only  performance  penalty  for  going  with  longer  signals  is  a  slight  increase  in  latency  at  the 
matched  filter,  and  the  need  for  filters  matched  at  every  scale.  In  interpreting  the  FOOBP  plot  of  figure  4, 
the  T  normalizing  the  horizontal  axis  is  the  reciprocal  of  the  transmission  rate  (in  bits/sec):  T  =  l/R. 


3  Fast  Rician  Fading  Channel  Model 

For  the  purposes  of  analysis  through  the  fading  channel,  the  signal  of  (11)  is  assumed  to  be  binary,  6^,/  €  ±1. 
The  model  for  the  fast  fading  Rician  channel  is  shown  in  figure  5  [9],  where 

Sd(t )  =  Re{v/25</(t)exp[7'2*/ct]} 

and  (/(f)  is  given  by  (11).  The  signal  sr(t)  is  the  fading  component,  given  by 

sr(f)  =  Re[£(t)<7(t  -  tD)  exp[j2ir(fc  +  fD)(t  -  fD)]  (13) 

where  tD  is  the  delay  of  the  fading  component,  fD  is  the  Doppler  shift  of  the  signal,  and  f(t)  is  a  complex, 
stationary,  zero-mean  Gaussian  random  process  with  autocorrelation  function 

Ri(r)  =  \E[C(mt  +  r)}  =  DPi(r). 

1)  is  a  constant  and  p<r(r)  is  normalized  so  that  p^(0)  =  1.  The  baseband  equivalent  fading  component  is 

sr (t)  =  t(t)g(t  -  tD)  exp[j2ir(f Dt  -  (fc  +  /d)7>)]  (I4) 

A',  +  1 

=  (;{t)  -  tD  -  iT'7’)exp\j2n{fDt  -  {fe  +  fD)tD)}-  (15) 

<7- 1  t£Z 

The  channel  also  introduces  the  AWGN  baseband  noise  process  n0(£)  with  correlation 

fin(r)  =  i£K(<)n0(t  +  r)]  =  N0S(t). 

The  baseband  equivalent  received  signal  is 

r(£)  =  sd{t)  +  sr(£)  +  n0(t). 

The  receiver  uses  filters  matched  to  the  signal  on  each  scale  and  employs  differential  detection.  A 
baseband  equivalent  receiver  is  shown  in  figure  6(a)  6(b).  The  receiver  for  scale  s  correlates  the  received 
signal  with  a  signal  p8{t  -  kT)  —  the  matched  waveform  starting  at  the  Arth  signal  interval  for  k  €  2s  Z  — 
and  integrates  over  [ kT ,  (k  +  qa)T}. 

The  output  of  the  matched  filter  on  scale  s  at  the  sample  instant  is 


z,((k  +  qB)T)  =  ztil  ( k )  =  c,((k  +  q,)T )  +  r,((*  +  qa)T)  +  n,((k  +  qeT) 

where  ca  is  the  direct  signal  part,  r,  is  the  fading  signal  part,  and  na  is  the  noise  part.  The  signal  part  can 
be  written  as 


ca((k  +  qa)T)  =  c,A 


(k) 


=  /“■ 


( k+q.)T 


ps{t  -  kT)sd{t)dt 
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which,  due  to  shift  and  scale  orthogonality,  may  be  written  as 

ct,i(k)=b.,ky/2S. 


The  noise  part  of  the  output, 

f(k+q.)T 

n8((k  +  q»)T)  =  n,,i  =  /  P.(*  -  kT)n0(t)dt 

JkT 


is  zero-mean  Gaussian  with  variance 

n.,i  =  ^E[n,(i 7,  +  kT)ns(q ,  +  *T)]  =  JV0,  (16) 

due  to  the  unit-energy  of  the  waveform  p,(£).  The  fading  component  output  of  the  matched  filter  is 

r(k+q.)T 

r,((k  +  q,)T)  =  /.,i(fc)  =  /  P*(t  -  kT)sr(t)dt. 

JkT 

For  differential  detection,  the  matched  filter  output  at  the  previous  sample  interval  is  also  used,  which  is 
z,((k  -  2‘  +  q,)T )  =  c8((k  -  2*'  +  q„)T)  +  r8((k  -  2*'  +  q,)T)  +  n8((k  -  2*'  +  q,)T) 


for  A-  €  2*’Z.  The  direct  signal  part  is 


r,  ((It  -  2*'  +  q.)T)  =  c,a{k)  =  f{  ^  p,(t  -  (k-  2  >')T)sd{t)dt  = 

J{k-2‘’)T 

Th»*  noise  part  n8 12  again  has  variance  The  covariance  between  n5> i  and  is 

1  r(Q.-r)T 

n  12  =  =  No  +  2®  T)dt  —  n^i  —  0. 

2  7o 

The  fading  component  from  the  previous  symbol  instant  is 

r,((Ar  -2a'  +  q,)T)  =  f.-2{k)  =  [  pa{t  -  (k  -  2s' )T)sr(t)dt. 

J (k— 2*'  )T 


(17) 


3.1  Probability  of  Error 

Although  the  symbols  at  other  scales  and  shifts  are  orthogonal  to  the  symbol  transmitted  at  scale  s  and 
interval  k  transmitted  in  the  direct  portion  of  the  received  signal,  the  random  multiplicative  factor  and  delay 
of  the  fading  component  sr(£)  makes  it  so  that  the  elements  of  sr(£)  are  not  necessarily  orthogonal  to  signal 
at  scale  s  or  to  each  other. The  probability  of  error  must  account  for  this  interference  for  each  possible  bit 
pattern  in  the  interfering  signals  and  for  each  alignment  of  the  matched  filter  relative  to  the  other  scales. 

Let  I,  k  be  the  set  of  bits  on  other  scales  that  overlap  with  the  signal  on  scale  s  and  starting  position  k 
and  let  /',  *('  !!,.*)  denote  the  probability  of  (bit)  error  for  the  Atth  position  symbol  at  scale  s,  given  all  the 
interfering  bits.  Note  that  it  is  necessary  to  express  the  probability  as  a  function  of  the  symbol  interval  k, 
since  amount  of  overlap  of  the  matched  filter  on  scale  a  with  other  symbols  depends  upon  the  position  of  the 
matched  filter  relative  to  the  symbols  on  the  other  scales.  The  conditioning  on  the  bits  may  be  eliminated 
by 


P.,k(e) 


1 

\I,k\ 


£P.,*(e|I ..*), 
x..» 


(18) 


wheie  the  sum  is  taken  over  all  possible  interfering  bits  and  \I\  denotes  the  cardinality  of  the  set  X.  The 
probability  of  error  for  symbols  at  scale  s  is  obtained  by  averaging  the  probability  of  errors  over  the  different 
matched  filter  alignments, 

P»(e)  =  ^37  P»,2*'i(e)- 

z  «=0 
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Finally,  the  overall  probability  of  error  is  a  combination  of  the  probabilities  on  each  scale,  weighted  by  the 
fraction  of  bits  in  the  overall  transmission  that  cure  sent  on  each  scale, 

p(e)=  £  ^7P’(e)'  (1Q) 

Conditioned  upon  the  interfering  bits,  the  fading  components  of  the  matched  filter  outputs  are  Gaussian 
and  the  vector 


/«,  i(*)  +n«,  1 

fs;2{k)  +  ti8,2 

is  Gaussian  with  zero  mean  and  covariance 


v(*|Z.,t)  = 


za,i(k)  -  c„,i(/r) 

~  c»Mk) 


K,(k\la,k)  =  |E[v"(A:|Z,tfc)v(*|Z.1*)]  = 


*..n(*)  K.,i2  (k) 

K.Mk)  K.,  2a(*) 


Assuming  that  no(<)  and  £(t)  are  statistically  independent,  then 


r,(*;l,l|Z,.k)  rg(A:;  l,2|Zs,jfc) 

r.(*;2,l|X.,*)  r,(*;2,2|Z,ifc)  J  ’ 

where  the  correlation  of  the  fading  components  is 

rt{k;n,m\llik)  =  ^ E[rt<n(k)r;>m(k)\I,tk }. 

In  the  appendix,  expressions  for  the  correlation  of  the  fading  components  are  developed. 

Let  0,  i t  he  the  difference  in  angle  between  za,i{k)  and  za^{k).  If  the  bits  at  interval  k  and  k  -  28'  axe 
the  same,  then  a  correct  decision  is  made  if  -n/2  <  9a  k  <  7t/2.  The  probability  of  this  event  not  occurring 

is  [9] 


Ka(k\la,k)  = 


ii  nati2 

n8,2i  n8  22 . 


+ 


P.,k(e\l»,k)  =  ~  F.,k(x/ 2), 

when*  b n  k  is  given  below.  If  the  bits  at  interval  k  and  k  —  2s  are  of  opposite  sign,  then  a  correct  decision  is 
made  if  rr/2  <  9a >  <  3tt/2  and  the  probability  of  error  is 


P*Ae\I,,k)  =  F9yk(n/ 2)  -  Fm(3tt/2)  =  FJt*(7r/2)  -  F,,k(- tt/2). 


The  function  /v*  is  defined  by  [10] 


where 


£  C,>  f _ M  sin(A^>aJt  -  j&) _ 

4tr  1C,,*  -  V,  ,*sin<  -  IT,,*  cos(A$gi*  -  rjj)cost 
*78,*  sin  ^  —  A,, t  cost/)  j 
•  -  cos^  +  As,*  sintA)  cost  J 


(20) 


C,,*  -  V',,*  sin  f  -  IT,  *  cos(A<i>,,*  -jp)cos t 
1  —  (t7s,*  cos  v  4-  A,,*  sin  ip)  cos  t 


^  8,*  —  y/Pl^TjcPi^Tji 


=  arg[c,)2(*)]  -  arg[cSii  (A:)] 

=  2^P'sk  +P2.».*)  C,,*  =  -(/>2,s,Jfc  7*1,8,*) 
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Be{K.Mk)} 

fKZi(k)K.,n{k) 


fc  — 


Im{if,,i2(fc)} 

'*..n(*)tf.,2a(*) 


Pl,s,k  = 


|c»,i(fc)|2  K2WJI 

2K,n(k)  P2'S ’  2Kt, 22(h) 


For  the  binary  differential  signaling  described  in  [9],  substantial  simplifications  result  in  these  equations. 
However,  for  multiscale  signaling,  such  simplification  is  not  possible  so  the  results  presented  below  ar 

obtained  by  numerical  integration  of  (20).  .  .  .  .  , 

To  compute  F.  k  it  is  necessary  to  know  the  correlation  matrix  K,(h\l.,k)  which  requires  knowing  the 

correlations  of  the' fading  components-  Expressions  for  the  fading  components  are  given  in  (79)  and  (80) 
in  the  appendix.  Due  to  the  complexity  of  the  expression  for  the  correlation  of  the  fading  component  and 
the  number  of  potential  interfering  bits  it  is  not  computationally  feasible  to  compute  the  correlation 

for  every  possible  interfering  bit  pattern.  An  approximation  is  therefore  made  to  the  probability  of  err 
bv  averaging  the  value  of  Fg<k  obtained  over  500  randomly  generated  bit  patterns.  (It  was i  determined  b> 
simulation  that  the  probability  of  error  calculated  by  averaging  F„tk  obtained  from  500  and  10000  random  y 

generated  bit  patterns  are  basically  the  same.  )  -  ,  ,  .  Tn 

Plots  of  probability  of  error  have  been  generated  to  demonstrate  the  performance  of  the  wavelets.  In  these 
plots,  the  correlation  function  is  for  a  land-mobile  fading  channel,  R^t)  -  DJ0{2ttBdt).  The  particular 
ra.se  of  L  =  0  t n  =  0  and  S/D  =  10  dB  is  shown  in  the  probability  of  error  plots.  As  a  summary,  there  are 
four  general  observations.  First,  for  both  small  and  large  BDT,  single  scale  wavelet  signalings  (signalings 
onlv  with  the  scaling  functions)  generally  perform  better  than  conventional  (single-scale  flat-topped  pulse, 
as  reported  in  [9]  )  signal.  Second,  there  is  very  little  change  in  the  probability  of  error  among  single  scale 
wavelet  signals  with  different  support.  Third,  the  penalty  for  using  multiscale  signaling  in  term  of  increase 
in  probability  of  error,  reduces  for  using  wavelets  with  longer  support.  And  fourth,  the  penalty  from  one 

scale  to  next  diminishes  with  the  increases  of  the  number  of  scales. 

Figure  7  illustrates  the  probability  of  error  for  BDT  =  0.01  for  conventional  and  wavelet  signaling^ 
In  figure  7(a),  wavelet  signalings  are  with  one  scale.  It  can  be  seen  that  there  is  approximately  a  1  dB 
improvement  of  wavelets  over  conventional  signaling.  There  is  very  little  change  in  the  probability  of  error 
among  D<  lh  D6  and  D8.  Figure  7(b)  illustrates  the  probability  of  error  for  four  scales  of  multiscale 
signalings*'  From  the  plots,  one  can  observe  that  D4,  De  and  D8  perform  just  as  well  as  conventional 
signaling.  However,  there  is  a  huge  penalty  for  D,  due  to  the  increase  in  correlation  among  signals  from 

different  scales.  .  .  c 

Figure  8  illustrates  the  probability  of  error  for  BpT  =  0.2.  The  basic  observation  is  the  same  for  figure 
8(a)  as  for  figure  7(a).  In  figure  8(b),  it  can  be  observed  that  the  penalty  is  larger  when  BpT  =  0.2  then  that 
is  when  II p  i  =  0  01  for  all  the  families  of  wavelet  except  for  Dt.  Actually,  D2  with  four  scales  is  performing 
better  in  channels  with  BDT  =  0.2  then  in  channel  BpT  =  0.01  while  all  other  wavelets  are  performing  not 

^Figure  9  illustrates  the  probability  of  error  for  BpT  =  0.5.  A  diminishing  increment  of  penalty  as 
the  number  of  scales  increases  is  observed.  In  figure  9(a),  while  D2  with  single  scale  has  around  1  d 
improvement  over  conventional  signaling  through  out  the  range  shown,  all  signals  with  one  scale  for  other 
families  of  wavelets  start  out  with  a  better  performance  when  EB/N0  is  low  and  become  having  not  as 
good  a  performance  then  conventional  signaling  as  EB/N0  getting  bigger.  The  crossing  begins  at  around 
Eb/.Yo  =  12  dB.  first  for  D8,  then  for  £>«,  and  then  for  D4  at  around  EB/N0  =  14  dB  The  above  order  also 
shows  the  increasing  performance  through  different  wavelet  families.  As  for  two  scale  signals,  figure  9(b) 
illustrates  that  there  is  a  penalty  for  all  wavelet  signalings.  D2  is  performing  as  conventional  signaling,  whi  e 
all  other  wavelets  do  not  perform  as  well.  D8  being  a  wavelet  having  a  longer  support  then  D,  suffers  less 
performance  penalty.  When  the  number  of  scale  increases  to  three,  figure  9(c)  shows  that  as  the  increment 
of  penalty  is  the  highest,  D2  becomes  the  worst  performer.  Now,  the  order  of  lower  probability  of  error 
becomes  conventional  signaling,  follows  by  D8,  D4,  Dt,  then  D2.  As  the  increment  of  penalty  diminishes, 
the  plots  of  probability  of  error  for  four  scales  signaling  in  figure  9(d)  is  nearly  the  same  as  that  in  figure 

9(c). 

From  figure  10,  it  is  observed  that  that  the  performance  variation  as  the  number  of  scales  increases  in 
smallest  for  the  D8  signaling. 
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3.2  Conclusions 

As  this  report  has  demonstrated  to  this  point,  wavelet  signals  can  be  used  for  bandwidth  efficient  communi¬ 
cation.  Most  of  the  benefits  are  obtained  when  single  scale  scaling  functions  are  used  for  transmission.  For 
bandwidth  consideration,  the  longer  the  support  of  the  wavelet,  the  more  bandwidth  efficient.  For  signals 
over  the  fast  Rician  fading  channel,  the  performance  does  not  seem  to  depend  significantly  upon  the  type  of 
wavelet  used,  with  shorter  wavelets  providing  only  slightly  more  improvement  than  longer  wavelets.  How¬ 
ever,  if  multiscale  wavelet  communication  is  contemplated,  there  is  no  change  in  bandwidth  efficiency  over 
single  scale  wavelet  communication.  But  for  the  fading  channel,  there  is  a  penalty  in  terms  of  increasing 
probability  of  error  for  multiscale  signaling.  This  penalty  decreases  as  wavelets  with  longer  support  are  used. 
At  the  same  time,  there  is  also  diminishes  with  the  increment  of  the  penalty  when  the  number  of  scales  is 
higher. 

Because  the  improvement  is  only  modest,  there  is  motivation  to  consider  signal  designs  which  are  optimal 
for  fading  channels.  Aspects  of  the  analysis  and  design  methodology  in  this  regard  axe  thus  taken  up  in  the 
next  Part. 


Figure  !:  Probabilitv  of  error  for  differential  detection  of  BPSK  for  different  fading  bandwidth  BpT  (from 

[Op- 
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^2T(0  =  Pl(0 


Figure  2:  (a)  Illustration  of  the  time  overlap  of  multiscale  transmission  for  Na  =  2  and  q  =  5.  (b)  A  matched 
filter  interval  for  scale  1  that  is  scale-aligned  with  scale  2.  (c)  A  non-scale-aligned  matched  filter  interval 


i>2  T 

T 

047- 


Figure  3:  Relative  placement  in  time  of  successive  symbols  on  multiple  scales. 
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Figure  4:  Fractional  out-of-band  power  (FOOBP)  for  QPSK,  MSK,  GMSK,  and  scaling  functions  with  4,  6, 
8,  and  10  coefficients. 


Figure  5:  Rician  fading  channel  model 
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Differential  dector  for 
Pi 


j 


Differential  dector  for 

PVS+ 1 

(a)  Multiscale  differential  detector 


Sampled  every  q8>  T 


k  such  that  k  €  2s  Z  and  k  <  q8'  <  (k  4-  1) 


(b)  Differential  detector  for  individual  ps 


Figure  6:  Differential  detection  for  multiscale  signaling 
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(a)  Comparison  of  probability  of  error  for  conventional  (b)  Comparison  of  probability  of  error  for  conventional 

signaling  and  single  scale  wavelet  signalings  and  wavelet  signalings  with  four  scales 


Figure  7:  Probability  of  error  comparisons  for  B^T  —  0.01 


4  A  Matched  Filter  Bound  for  Wavelet  Transmission  Through  a 
Fading  Channel 

In  this  section,  an  analysis  is  provided  of  the  performance  of  multiscale  transmission  (wavelet  modulation) 
using  the  formalism  of  the  matched  filter  bound.  The  matched  filter  bound  [11,  12]  is  a  lower  bound  on 
performance,  typically  employed  in  a  fading  channel,  which  characterizes  the  detection  capability  of  a  single 
communication  pulse  transmitted  in  isolation,  so  that  potential  effects  of  ISI  are  neglected.  Even  though  the 
match  filter  bound  (MFB)  is  an  optimal  lower  performance  bound  over  communication  channel  that  may 
not  be  practically  realizable,  it  is  a  common  performance  measurement  for  communication  systems  [11,  12]. 
There  are  reasons  for  saying  MFB  is  an  optimal  performance  bound.  First,  it  is  assumed  that  all  the 
information  about  the  channel  is  known,  which  is  not  feasible  for  a  real  fading  channel.  Second,  transmitted 
pulses  are  separated  sufficiently  so  that  no  intersymbol  interference  (ISI)  occurs.  In  other  words,  the  error 
rate  is  determined  by  assuming  only  one  data  pulse  is  transmitted. 

In  the  present  case,  consideration  is  given  to  a  set  of  pulses  transmitted  at  different  scales  over  a  common 
time  interval,  but  otherwise  isolated  from  other  pulses.  The  MFB  for  multiscale  wavelet  pulse  is  derived 
in  section  4.2.  In  section  4.3,  the  MFB  is  presented  for  examples  from  Daubechies’  wavelet  families  and 
compare  with  that  of  raised  cosine.  Section  4.4  presents  conclusions  and  discussions. 

4.1  Channel  Model  and  Wavelet  Transmission 

As  before,  the  signal  is  modeled  as 


N,-l 

9(0 =  y  ]  baPa(0m  (21) 

a-  1 

For  the  MFB  analysis  of  this  section,  only  one  wavelet  function  pulse  is  present  on  each  of  the  different 
scales. 

A  baseband  equivalent  Rayleigh  fading  channel  consisting  of  p  paths  can  be  represented  as 

p 

c(t)  =  Y^aizi(t)6(t -Ti),  (22) 

*=1 
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(a)  Comparison  of  probability  of  error  for  conventional  (b)  Comparison  of  probability  of  error  for  conventional 

signaling  and  single  scale  wavelet  signalings  and  wavelet  signalings  with  four  scales 


Figure  8:  Probability  of  error  comparisons  for  BoT  =  0.2 


where  z,{t)  is  a  slowly  time- varying,  zero-mean,  unit  variance  complex  Gaussian  random  process,  tj  is  the 
delay  of  the  zth  path,  and  a,  is  the  root-mean -square  value  of  the  magnitude  of  ith  path  [11].  Assuming  that 
z,(t)  does  not  change  within  the  duration  of  a  single  pulse,  z<  can  be  written  in  place  of  z,(<).  The  baseband 
signaling  pulse  g(t)  is  sent  through  the  channel  (22)  and  corrupted  by  additive  white  Gaussian  noise.  At  the 
receiver,  the  received  signal  is 

r(t)  =  ]T]  aiZi9(t  -  Ti)  +  n{t).  (23) 


4.2  Derivation  of  Matched  Filter  Bound 

To  determine  the  data  sent  with  multiscale  signaling  pulse,  pass  the  received  signal ^(23)  ^through  a  bank  of 
matched  filters  as  shown  in  figure  11.  For  each  scale  £,  r(t)  is  matched  to  hPc(t)  =  £f=i  oii  zi  p^[T^  —  {t  —  r,)], 
then  the  signal  is  sampled  at  t  =  to  give  us  the  received  signal  rp<1  where  Tq  is  the  support  of  pc(t). 
Now.  let  us  derive  the  error  rate  P(e<;)  for  the  data  transmitted  in  scale  £.  The  value  of  rp<  can  be 
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(a)  Comparison  of  probability  of  error  for  conventional 
signaling  and  single  scale  wavelet  signalings 


(c)  Comparison  of  probability  of  error  for  conventional 
and  wavelet  signalings  with  three  scales 


(b)  Comparison  of  probability  of  error  for  conventional 
and  wavelet  signalings  with  two  scales 


(d)  Comparison  of  probability  of  error  for  conventional 
and  wavelet  signalings  with  four  scales 


Figure  9:  Probability  of  error  comparisons  for  BqT  =  0.5 
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expressed  as 


'p< 


=  fT<  r(t)hP((T(-t)dt 
Jo 

=  [  E  ocjZjg(t  -  n)  +  n(t)]E  ajZjP^it  -  7j)]<ft 

■'°  t=i  i=1 

=  fT<  e  Q.^i  e'  -  r^dt 

Jo  <=1  <r= 1  j=l 


T  P 

+  [  '  nit^a-z’pKt-rj)}* 

Jo  j= i 

JV.+l  P  P  /•?< 

=  e  { EEnriQ72,2^<T  /  — t»)pc^ — 

+  EQj2j  [  n(f)Pc(f  -  Ti)dt 

j=  i  Jo 

=  e' b° - r>} +EQi2i l c n(t)p<(t - T')d* 

<r  =  l  i=  1  j  =  l  J-1 

=  E  KZHMaiz  +  E QJ2J  /  n(f)Pc(*  - 


(24) 


£7  --  1  J— 1 

where  Z  -  [*„••■  ,;P]T  and  MctC  =  m.Cij  =  0,0$  JW*  "  Tf)-  Assume  6<  =  L  Let  6  = 

.1....  ,6jv..i],  and  the  set  of  all  possible  b  be  3.  Note  that  \B\  =  2N‘.  Conditioning  on  6,  it  ■» 

found  that 


is 


N.  +  l 


£J(rp<:  |6)  =  ^2  b(rZH Ma(Z  -  EPc , 


a- 1 


and 


So,  the  pdf  of  rK  is 


For  Ev ^  >  0 


4  =  yz"m<z- 


,,  I,  N  1  r  (rP<  fipc)’  1 

/(rpJfe)  =  -7^—  exP[ - 2^ - ] 


\/2ncr 


(25) 


(26) 


(27) 


PC 


PC 


P(ec\b)  =  J  f(rPc\b)drP( 

1  f°  r  (rPC  “  EPc)2,  . 

=  -7= -  /  exp[ - - ]®p< 

V27T crpc  y~oo 

Ep< 

_i__  r~p 

n/2tt  7—00 


(28) 


—  i—  j 

e  2  dx 


yPC 
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For  EP(.  <  0,  by  the  same  argument, 


P(ec\b)  =  l-Q 


Pc 


(29) 


Also,  if  =  —  1, 


P(e<\b)  = 


_M£J 


FPc  <0 


1  -  Q  [~^j  otherwise 
The  average  bit  error  probability  given  b  can  be  found  by  averaging  P(e^\b)  over  all  EP(: 


P«ve(\b  =  £°  P(e<\b)p(EP()dE,e. 


Conditioning  on  b  is  eliminated  by 


Pave^  |g|  ^>ave<\b' 

bee 


(30) 


(31) 


(32) 


Finally,  as  one  bit  is  sent  through  each  of  the  scales,  the  average  error  rate  over  the  matched  filter  bank  is 

tf.  +  l 


Pave  ~  N,  +  1  ^  Pave(' 


(33) 


4.3  Examples 

In  this  section,  the  results  of  the  previous  section  are  used  to  derive  the  MFBs  of  some  multiscale  wavelet 
signalings  for  some  practical  fading  multipath  models.  The  channel  models  used  are  typical  urban  (TU) 
and  hilly-terrain  (HT),  which  axe  used  as  benchmarks  in  the  Pan-European  Digital  Cellular  (GSM)  standard 
[11].  The  wavelet  functions  used  here  are  from  Daubechies*  wavelet  families,  in  particular,  D4,  Dq  and  Dg 
art'  list'd,  where  for  Dn,  n  is  the  number  of  coefficients  used  to  generate  the  wavelet  functions  [1],  These 
MFSs  will  be  compared  with  traditional  square-root  raised  cosine  signaling  with  excess  bandwidth  /?  =  0.35 
as  specified  by  IS-54  [11].  The  delay  and  amplitude  profiles  of  the  6-path  Rayleigh  fading  channel  models 
are  listed  in  Table  1.  Two  different  symbol  intervals,  one  for  IS-54  and  one  for  GSM,  are  studied. 

For  IS-54,  the  symbol  interval  used  is  40/xs.  Figure  12  illustrates  the  MFB  for  TU  channel.  From  the 
plots,  it  can  be  seen  that  all  MFBs  are  nearly  the  same  with  wavelet  signalings  performing  a  little  better 
than  raised  cosine.  Among  different  wravelet  signalings,  D4  has  the  best  MFB,  while  Dq  has  a  MFB  which 
is  the  same  as  that  of  the  raised  cosine.  In  considering  the  effect  of  different  number  of  signaling  scales,  it  is 
seen  that  the  MFB  for  tw*o  scales  of  D4  is  similar  to  that  of  with  one  scale,  wffiile  for  Dq  and  Dg,  two  scale 
signaling  gives  marginally  better  MFB  over  one  scale  signaling.  Figure  13  illustrates  the  MFB  for  the  HT 
channel.  Again,  the  MFB  for  single  scale  wravelet  signalings  are  similar  to  those  of  raised  cosine.  However, 
for  all  D4  to  Dg,  twx>  scale  signaling  have  MFBs  inferior  to  that  of  single  scale  signaling. 

For  GSM,  the  symbol  interval  used  is  3.8/is.  As  the  symbol  interval  decreases,  it  can  be  seen  by  comparing 
figures  14  and  15  to  figures  12  and  13  that  there  is  a  general  improvement  in  probability  of  error.  Figure  14 

illustrates  the  MFB  for  TU  channel.  It  can  be  seen  that  for  one  scale  signaling,  D4  has  a  2 - 3  dB  gain 

over  raised  cosine.  This  gain  decreases  wffien  longer  w’avelet  is  used.  However,  as  the  MFBs  for  HT  channel 
illustrated  in  figure  15,  all  single  scale  wavelet  pulse  perform  the  same  as  raised  cosine.  For  both  TU  and 
HT  channels  w'ith  GSM  symbol  interval,  increasing  the  number  of  scales  for  wavelet  pulse  from  one  to  two 
worsens  the  MFB. 

4.4  Conclusions  and  Discussions  Regarding  the  MFB 

In  this  section,  the  matched  filter  bound  (MFB)  of  time-discrete  multipath  Rayleigh  fading  channels  derived 
in  [13,  11]  was  extended  to  multiscale  wavelet  signaling  communication.  From  the  examples  presented  in 
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ale  wavelet  pulses  have  better  MFBs  compare  to  rarsed 
section  4„«  c-  U— ~  ~  Si"6e;  FB  roeva.^dl^cm— — 

C^ltowever,one  tlm^dnhe^jorma^j^^g  ^^^ngwith 

symbol  intervals  are  use  , 

^  ...  .  profiie  of  GSM  ChannelstU]- 

Table  1:  Delay  and  Amphtu 

HT  model 


Path  # 
1 
2 

3 

4 

5 

6 


TU  model 

Delay  (ps)  Po^’er 

0.000  I-000 


0.813 

1.626 

2.439 

3.252 

4.056 


1.000 

0.669 

0.448 

0.300 

0.200 

0.134 


HT  model 

Delay(ps)  Power 

o.ooo  i-000 


0.813 

1.626 

15.447 

16.260 

17.073 


1.000 

0.251 

0.060 

2.258 

0.177 

0.122 


5  Amplitude  vs.  Phase  Modulation^ 

method  beOT 


6  Spectral  Adaptivity  diffae„,  time  and *■»-£££££; 

:rr^ 

*“*  ^  Irf  the  nofse  spectrum  ot  the  chance!  |IA  1  Noru  „  Utah  State 

:§sssSu^-— 8 

be  found  [18]. 
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Average  Probability  of  Error 


Figure  10:  Illustration  of  probability  of  error  penalty  for  multiscale  wavelet  signalings 
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Figure  11:  Matched  Filter  Bank  for  multiscale  wavelet  transmission 
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Average  Probability  of  Error  Average  Probability 


Et/Vdb> 


Eb'No<db> 


Figure  12:  Typical  urban  with  symbol  interval  40 (i  s 
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Figure  13:  Hilly-terrain  with  symbol  interval  40 [i  s 


Part  II 

Fading  Channels 

7  Background  and  literature  review 

At  present,  wireless  communication,  mainly  cellular,  systems  are  only  able  to  provide  low  speed  data  com¬ 
munication  at  error  rates  that  are  far  from  acceptable  in  wireless  network  connections.  This  is  because  these 
systems  are  based  on  technology  distant  from  theoretical  limits  [19].  One  of  the  main  reasons  for  this  draw 
back  in  technology  is  the  hurdle  of  fading  multipath  phenomenon  [20].  For  example,  the  cellular  system, 
GSM,  we  are  using  today,  are  designed  more  to  circumvent  the  problems  imposed  by  the  fading  multipath 
phenomenon  than  actually  solving  them.  To  avoid  a  too  rapid  fluctuation  of  the  received  signal  power  (caused 
by  fading),  specific  signaling  rate  is  chosen  (according  to  the  rate  of  fading)  so  that  the  received  signal  can  be 
regarded  approximately  constant  over  a  data  symbol  interval  (slow  fading  assumption).  At  the  same  time, 
the  signaling  rate  can  not  be  too  high  as  that  will  cause  severe  problems  with  Inter  Symbol  Interference,  ISI 
(frequency  selective  due  to  multipath  phenomenon)  [21].  However,  the  good  thing  is  there  is  an  indication 
that  high-speed  reliable  communication  is  possible  if  the  user  is  willing  to  make  the  most  of  what  the  media 
would  give  him  [22].  Because  of  this  belief,  considerable  research  effort  has  been  devoted  in  the  area  of 
multipath  fading  channel  communication  for  the  past  five  decades  [23,  24,  25,  26,  27,  28,  29,  30,  31].  This  is 
because  fading  multipath  phenomenon  is  unavoidable  for  high-speed  wireless  communication. 

Tsually  fading  multipath  channels  are  merely  called  fading  channels  [14].  However,  there  are  actually 
two  basic  degradations:  i.e.,  multipaths  and  deep  fades.  The  multipaths,  which  results  from  reflections, 
that  cause  transmitted  signals  being  superimposed  at  the  receiver.  When  the  arrival  time  of  the  different 
rays  are  of  the  same  order  of  magnitude  as  the  duration  of  the  transmitted  signals,  successive  signals  are 
smeared  together,  and  cause  ISI.  The  span  of  the  excess  delay  is  directly  related  to  the  distortion  rate 
of  the  channel  to  the  signal  according  to  the  difference  in  frequency.  Frequency  selectiveness  is  the  term 
we  use  when  the  distortion  rate  changes  a  lot  for  a  small  variation  in  frequency.  When  the  arrival  time 
difference  is  comparable  to  the  period  of  the  carrier  frequency,  another  effect  results,  i.e.,  the  deep  fades. 
The  superposition  of  many  waves  of  different  phases  give  a  spatial  interference  and  cause  extremely  low 
signal  power.  Deep  fades  can  also  be  caused  by  relative  motion  between  the  transmitter  and  the  receiver, 
which  also  results  in  superposition  interference  that  is  called  the  Doppler  Effect  [32].  Of  course,  deep  fades 
can  also  be  caused  by  the  continuous  physical  changes  of  the  transmission  media.  When  the  interference 
and  fluctuation  are  a  lot  more  rapid  compared  with  the  signaling  rate,  fast  fading  occurs. 

A  lot  of  time,  fading  and  multipath  phenomenon  are  coming  hand  in  hand.  However,  there  are  situations 
like  stationary  communications,  i.e.,  both  the  transmitter  and  the  receiver  are  not  moving,  where  the  occur¬ 
rence  of  multipath  transmission  is  the  main  cause  of  degradation  of  communication  quality;  and  there  are 
also  situations  like  communication  betwreen  high  speed  vehicles  at  an  environment  with  no  major  reflection 
surfaces,  where  the  main  cause  of  degradation  is  fading.  To  make  the  multipath  fading  channel  problem  to 
be  more  tractable,  it  is  helpful  to  consider  only  the  multipath  aspect  or  the  fading  aspect  alone.  In  this  part 
of  the  report,  only  the  fading  aspect  is  studied. 

7.1  General  Fading  Channel  Model 

Fading  channels  w*ere  first  modeled  as  large  number  of  “scatterers”  located  at  random  points  within  the  prop¬ 
agation  path  in  1950s  and  1960s  [25,  27].  This  idea  was  primarily  applied  to  over-the-horizon  (troposcatter) 
communic  ations  covering  a  wide  range  of  frequency  bands.  To  reach  beyond  the  horizon,  clouds  of  particles 
in  the  troposphere  were  used  as  reflectors  for  the  radio  waves.  The  differences  in  path  length  between  the 
large  number  of  scattered  weaves  give  rise  to  Rayleigh  fading  if  there  is  no  dominant  direct  component  [33]. 

The  ground  w’ork  of  modeling  fading  channels  mathematically  was  laid  by  Price  [23,  24].  For  the  following 
half  a  century,  there  w'ere  tremendous  amount  of  contributions  to  this  field.  In  a  recent  survey  paper  written 
by  Biglieri,  Proakis,  and  Shamai,  more  than  five  hundred  entries  of  citations  reporting  the  state-of-the-art 
achievement  in  the  research  of  fading  channel  wTere  included  in  their  reference  [31].  By  considering  the 
absence  of  Prices'  papers  in  their  citation,  how  many  more  others  were  “missing”!  Among  those  that  were 
not  included  in  their  citations,  there  was  a  paper  written  by  Bello  who  introduced  a  simple  way  to  model 
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the  fading  phenomenon  with  the  notion  of  a  wide-sense  stationary  uncorrelated  scattering  (WSSUS)  [26]. 
Uncorrelated  scattering  means  the  attenuation  and  phase  shift  of  the  channel  associated  with  different  path 
delay  are  uncorrelated. 

A  widely  accepted  mathematical  modeling  of  fading  channel  nowaday  is  presented  in  Proakis’  book  [30]. 
A  fading  channel  is  be  viewed  as  a  continuous  time- variant  filter  with  a  baseband  equivalent  impulse  response, 


(34) 


where  c(.; .)  can  be  complex,  r  is  the  delay  in  response  and  t  is  the  time  variable.  If  a  signal  s(t)  is  sent 
through  this  channel  with  Additive  Gaussian  White  Noise  (AGWN),  we  will  have  at  the  receiver  side, 


c(r;  t)s(t  -  r)dr  +  n(t), 


(35) 


where  r(t)  is  the  received  signal  and  n(f)  is  the  AGWN. 

Assuming  wide-sense  stationary,  the  autocorrelation  function  of  c(r;t)  can  be  defined  as 

<t>c{T\,T2\  Af)  =  ^£[c*(r1;t)c(r2;t  + At)].  (36) 

With  uncorrelated  scattering, 

<Pc(ti;  At)S{ri  -  To)  =  ^ E[cm (ti ;  t)c(r 2 ;  t  -t-  At)] .  (37) 

Setting  At  =  0,  the  expression  <^c(Ar;0)  in  (37)  is  usually  called  the  multipath  intensity  profile  or  the  delay 
power  spectrum  of  the  fading  channel.  Further  define  Tm  to  be  the  length  of  the  span  of  <£c(Ar;0)  that  is 
essentially  nonzero,  and  call  it  the  multipath  spread  of  the  fading  channel. 

Now.  if  we  take  the  Fourier  transform  of  the  channel  impulse  response, 

C{f;t)  =  j  c{T-,t)e-j2KfrdT.  (38) 

The  corresponding  autocorrelation  function,  with  WSSUS,  is 


OcU 1 ,  h\  At)  =  |  £[<?-(/, ;  t)C(fo\ t  +  At)] 

=  \jj  E[c*(n;t)c(r2:t  +  Af)]e-''2,r(/'r'-*r:')dr1dr2 
=  j  J <Mn;  At)<S(n  -  r2)eJ-’r(/'  ^-'^dndro 

=  J  <t>c(n-,&t)ej'2*u'-f')T'dTl 
=  J  <A0(r1;At)e-j2,rA/r‘dr1 

=  ^c(A/;Af) 


(39) 


Setting  At  —  (),  one  obtains  <pc(Af\0)\  it  is  the  autocorrelation  function  in  the  frequency  variable.  Define 
the  length  of  the  span  of  $c(A/;0)  that  is  essentially  nonzero  as  (A f)c-  Usually  (A f)c  is  denoted  as  the 
coherence  bandwidth,  which  is  a  measure  of  the  frequency  coherence  of  the  fading  channel.  To  illustrate 
this  idea,  assume  a  signal  with  a  bandwidth  less  than  (A/)c  is  transmitted  through  this  fading  channel. 
As  all  frequency  components  will  be  affected  essentially  similarly  by  the  channel,  this  channel  is  said  to 
be  frequenry-nonselective.  However,  if  another  signal  with  a  bandwidth  larger  than  (A f)c  is  transmitted 
through  the  same  fading  channel,  different  frequency  components  will  be  affected  differently  by  the  channel. 
At  this  time,  the  channel  is  said  to  be  frequency-selective.  In  this  case,  the  signal  can  be  severely  distorted 
by  the  channel.  Another  useful  result  due  to  the  fact  that  0c(Ar;O)  and  <fc(A/;0)  are  Fourier  transform 
pairs, 


(40) 
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To  study  the  time  variation  of  the  channel,  return  to  ^(A/;  At)  of  (39).  Setting  A /  —  0,  one  essentially 
obtains  the  autocorrelation  function  of  the  channel  in  time  averaging  through  all  the  different  frequency 
components.  Now,  the  length  of  the  span  of  <j>c{ 0;  At)  that  is  essentially  nonzero  is  defined  as  (At)c,  which 
is  called  the  coherence  time  of  the  fading  channel.  So,  for  a  digital  communication  system,  if  the  signaling 
period  is  shorter  than  (At)c,  it  may  be  assumed  that  the  channel  is  essentially  constant  for  individual  signals 
and  the  situation  of  a  slow  fading  occurs.  On  the  other  hand,  if  the  signaling  period  is  longer  than  (A f)c, 
it  may  no  longer  be  assumed  that  the  channel  is  constant  for  individual  signal,  and  fast  fading  situation 
occurs. 

In  practice,  (At)c  is  very  difficult  to  be  measured  directly,  so,  a  related  function  is  employed.  This  new 
function  is  called  the  scattering  function  of  the  fading  channel,  and  it  is  defined  as 

S(r;A)  =  J  J  <t>c(&f;At)e~2*TAfe-2*XAt  dAfdAt,  (41) 

where  r  is  the  time  delay,  and  A  is  the  Doppler  frequency  as  the  channel  time  variation  is  usually  caused 

by  physical  motions  of  the  media,  the  transmitter,  and  the  receiver.  The  scattering  function  of  the  channel 

provides  us  with  a  measure  of  the  average  power  output  of  the  channel.  When  r  =  0,  the  length  of  the  span 
of  5(0;  A)  that  is  essentially  nonzero  is  called  the  Doppler  spread  B<i  of  the  fading  channel.  Besides,  the 
following  relationship  exists  between  the  channel  coherence  time  and  the  Doppler  spread, 

(At )c  »  (42) 

Time  variation  autocorrelation  function,  0c(O;At),  can  be  obtained  from  5c(A)  s  5(0;  A),  called  the 
Doppler  power  spectrum  of  the  fading  channel.  Some  of  common  Spectra  and  autocorrelation  functions  are 
listed  in  Table  2  [34,  35,  36,  37,  9]. 


Table  2:  Spectra  and  autocorrelation  functions  of  fading  process 


Denotation 

Spectrum  5c (A) 

Autocorrelation  Function  Rc(j) 

1 .  Rectangular 

jk 

tv- sin  ‘InBnr 
'2-jvBdt 

2.  Gaussian 

Kexp{f£)/^BD 

K  exp[-(irBDT)2} 

3.  Land  Mobile 

K 

7r(A3-BD3)i/2 

KJo(2kBdt) 

4.  First-order  Butterworth 

K 

TbmTTTTbd7! 

K  exp[-2n Bd\t\] 

Note:  J0(-)  is  the  zero-order  Bessel  function  of  the  first  kind  and  K  is  a  constant. 

7.1.1  Categories  of  degradation 

Different  fading  phenomenon  introduces  different  degradations  to  communication  systems.  Those  degrada¬ 
tions  are  summarized  in  Table  3  [14]. 

7.1.2  Statistical  Models 

Let  us  revisit  (35),  i.e., 

r(£)  =  f  c(r;  t)s(t  -  r)dr  4-  n(t). 
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Table  3:  Categories  of  degradation  according  to  different  fading  phenomenon 


Frequency-nonselective 

Frequency-selective 

Slow  fading 

Fast  fading 

loss  in  SNR 

ISI  distortion, 
pulse  mutilation, 
irreducible  BER 

low  Doppler, 
loss  in  SNR 

high  Doppler, 
PLL  failure, 
irreducible  BER 

By  ParsevaTs  relation, 


r(t)  =  [  C(/;  t)S(f)ej2*ftdf  +  n(t).  (43) 

Assume  there  is  a  frequency  non-selectiveness  (that  C(f\t)  =  C(0;£)),  and 

r(t)  =  I  C(f-,t)S(f)e*’ft4f  +  n(t) 

=  <7(0;  t)  J  SU)e?3wft4f  +  n(t)  (44) 

=  C(0;t)s(f)  4-  n(t). 

The  transfer  function  C(0;f)  for  a  frequency-nonselective  fading  channel  may  be  expressed  as 

<7(0;t)  =  a(t)=a(*)e_J<W  (45) 

where  the  random  process  a(t)  represents  the  envelope  and  random  process  0(t)  represents  the  phase  of  the 
the  transfer  function.  In  this  case,  the  channel  is  also  called  a  multiplicative  fading  channel.  Please  note 
that  a  multiplicative  channel  can  be  both  fast  and  slow  fading.  For  different  fading  channels,  0(t)  is  always 
modeled  as  uniformly  distributed  over  the  interval  (-jr,7r).  Conversely,  for  any  fixed  value  of  t,  a(t)  can 
acquire  different  statistical  distributions,  from  which  different  fading  channel  is  named  after.  Among  common 
distributions  for  u(t)  are  the  Rayleigh  distribution,  the  Rician  distribution,  and  the  Nakagami  distribution. 

7.2  Literature  review 

All  researches  in  engineering  have  two  main  goals:  1.  understanding  and  2.  improvement.  Research  in  fading 
channel  is  not  an  exception,  two  characteristics  are  clear  in  till  of  the  works:  1.  finding  the  ultimate  performing 
bounds  for  fading  channel  (understanding),  and  2.  finding  ways  to  combat  fading  effects  (improvement).  An 
interesting  observation  is  at  the  early  days  of  the  study  of  fading  channel,  most  of  the  effort  was  devoted  to 
the  improvement  side  of  research,  only  until  later,  more  attention  was  attributed  to  understand  the  matter. 

The  history  of  fading  channel  studying  can  be  roughly  divided  into  three  periods.  The  first  period  spans 
from  the  early  1950’s  to  late  1960’s.  The  driving  force  behind  the  research  about  fading  channel  in  this 
period  was  mainly  the  development  of  long-distance  troposcatter  communication  [23,  24,  38,  33].  During 
this  period,  a  tremendous  effort  was  invested  to  collect  real  data  to  illustrate  the  fading  channel,  and  some 
practical  mathematical  models  were  built  [23,  26].  The  single  most  important  idea  for  combating  fading 
phenomenon,  i.e.,  “diversity”,  was  also  formed  in  this  period  [39,  40,  41,  42,  43,  44,  45,  46].  Even  a  lot 
of  work  was  performed  in  the  area  of  detection,  the  resulting  systems  were  not  reliable  enough  for  some 
computer  communication  applications.  This  was  the  main  reason  of  the  introduction  of  coding  theory  to 
signaling  on  the  fading  channel  [28]. 

The  second  period  spans  from  early  1970’s  to  mid  1980’s.  During  this  period,  there  was  a  drop  of  the 
interest  in  fading  channel  communications,  maybe  because  of  the  high  demand  in  the  research  of  wire  and 
optical  communications  where  fading  w-as  not  tin  important  issue.  Even  so,  there  was  still  continuous  effort 
in  the  Information  and  Coding  Theory  aspect  of  fading  channel  researches  [47,  48,  49,  50]. 
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The  third  period  spans  from  the  mid  1980’s  to  present.  During  this  period,  there  has  been  a  huge  increase 
of  scientific  activities  in  the  area  of  fading  channel  research.  The  main  driving  force  is  the  tremendous  increase 
in  the  demand  for  personal  mobile  wireless  communication,  of  which  the  following  references  typify  aspects 
of  implementations.  Many  new  implementation  schemes  have  been  derived  and  reported.  These  approaches 
can  be  further  subdivided  into  three  sub-classes,  i.e.,  diversity,  detection  and  estimation,  and  combined 
methods. 

1.  Diversity:  the  application  of  the  idea  of  diversity  exists  virtually  in  every  branch  of  research  of  fading, 
either  implicitly  or  explicitly.  A  basic  motivation  for  diversity  is  that,  under  a  fairly  general  condition,  a 
channel  affected  by  fading  can  be  turned  into  an  AWGN  channel  by  increasing  the  number  of  diversity 
branches  [46,  51].  There  are  mainly  two  types  to  diversities,  i.e.,  time  and  space  diversities.  Time 
diversity  is  usually  achieved  by  using  some  kind  of  rearrangement  of  the  signaling  [52],  or  by  means  of 
coding  [53,  51,  45,  54].  Space  diversity  can  be  achieved  by  using  multiple  antenna  at  the  transmitter 
and/or  the  receiver  side  [55,  56,  57]. 

2.  Estimation  and  Detection:  new  twists  to  adopt  old  techniques  to  fading  situations  [58,  59,  60].  Ap¬ 
plication  of  diversity  on  the  design  of  detector  and  estimator  [61,  46].  And  a  lot  of  reports  in  the 
application  of  equalizers,  besides,  for  a  lot  of  time  with  robust  adaptive  algorithm  for  fast  fading 
[62,  63,  64,  65,  66,  67,  68,  69]. 

3.  Combined  methods:  Information  theory  and  space  diversity  [70].  Adaptive  CDMA  signaling  [71].  And 
the  application  of  channel  coding  [72,  73]. 

8  Analysis  of  Multiplicative  Fading  Channel  Model 

This  section  considers  the  asymptotic  behavior  of  a  fading  channel  as  the  rate  of  fade  increases  compared 
to  the  support  of  signaling  function  when  a  correlation  filter  is  employed.  Because  the  study  deals  with  the 
behavior  of  the  channel  as  it  is  fading  very  fast  compare  to  the  length  of  the  signal,  it  can  be  assumed  that 
the  support  of  the  signal  is  divided  into  multiple  chips  and  the  signal  is  faded  independently  in  each  of  those 
chips.  This  will  greatly  simplify  the  development  of  the  arguments,  and  as  the  number  of  chips  is  large,  it 
should  reflect  some  interesting  channel  models. 

8.1  System  Model 

I  lje  system  model  is  presented  in  figure  16.  The  binary  information  a  =  {1,-1}  is  carried  by  the  signal 
g(t)  ~  us (f),  where  supp(s(t))  =  [0,T].  Then  g(t)  is  sent  through  a  multiplicative  fading  channel  with  a 
multiplving  factor  of  a (#),  with  a  being  positive  and  Rayleigh  distributed.  Also,  the  signal  is  corrupted  by 

No 

additive  Gaussian  white  noise  n(f)  with  variance  a  =  — .  At  the  receiver  side,  the  received  signal  is  first 

multiplied  by  .s(f),  then  integrated  through  the  support  of  s(t)  and  sampled  at  T  to  give  r.  If  r  >  0  a  is  1, 
else  a  is  -  1.  The  goal  is  to  evaluate  the  error  rate  versus  the  transmission  signed- to-noise  ratio.  To  simplify 

the  calculation,  it  is  assumed  that  s(f)  is  an  unit  energy,  i.e.,  J  s2(t)dt  =  1.  So,  the  detected  signal  r  can 

be  expressed  as 


(46) 


where  n  is  the  noise  component. 

For  fast  fading  situation,  divide  the  time  interval  of  the  support  of  g{t)  into  multiple  subintervals,  and 
assume  the  channel  varies  independently  in  each  subintervals.  Within  any  of  these  subintervals,  it  is  assumed 
that  the  channel  fades  slowly,  i.e.,  the  multiplicative  factor  Qj(t)  is  constant  in  each  subinterval,  so  at-  may 
be  employed  in  place  of  The  above  idea  is  illustrated  in  figure  17.  With  the  above  assumption,  (46) 

can  be  expressed  as 


r  =  r  4-  n 
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where 


=  X>t7?i 


and 


r*»+ 1 

m~L 


s2(t)dt 


8.2  Analysis  Approach  I 

Consider  that  the  interval  [0,T]  is  divided  into  N  chips.  Then  the  following  distribution  functions  for  the 
received  r  hold: 


/{4}(rlQ0i  Qi  i  •  •  •  ,  gn-i  ,  o  =  l)  =  — ^==e  (r 


and 


/{/»}(rlQo,Qi,  •  •  •  ,oA'-i,a  =  -1)  = . y=e  ^r+r'>^N°. 


n/ttA'o 

The  probability  of  error  given  q0,  Qj,  •  •  •  ,ccn-\  and  a  =  1  is 

f° 

p(e|Q0,Qi,---  ,ajv-i,a  =  1)  =  /  /{4}(r|a0,Qi, •  •  •  ,a/v-i,a  =  l)dr 
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i  r°° 

=  -L  f  e-‘l>dx 

n/27T 

and  similarly,  the  probability  of  error  given  a0,ai,  •  •  ■  7qs-i  and  a  =  -1  is 

fOO 

p{e\a0,aw--  ,Qn-i ,a-  -1)  =  /  f{A}(r\ot0,  Qi ,  *  *  *  ,Qyv-i,a  =  -1  )dr 

Jo 

i  r° c 

=  —j=  [  e~z3^'2dx. 

y/2 *  Jl 

If 11  =  I  and  a  =  —  1  are  equally  likely,  then  the  probability  of  error  given  Qo,Qi,-  •  •  ,qn-\  is 

, o/v-i)  =  ^p(e|Qo,Qi, •  •  •  ,Qn_i,a=  1)  +  ip(e|a0,ai,- •  •  ,aN-iya  =  -1) 

i  roo 

=  —f=  [  e~z3/1dx. 
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I  In1  probability  of  error  of  this  particular  detection  scheme  may  be  obtained  by  averaging  over  all  possible 
combinations  of  (a0,aj,  •  •  •  ,a/v_i),  i.e., 


{’ ~  J{  }  /{'t}(Qo.ai,- ••  ,oa-i)  |-^==y^  e  l3/2(ir|  d(Q0,Qi,--- ,qn_, 


)• 


Because  {o}s  are  independent,  it  follows  that 
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Let  Wt  =  Y'  — Qi ■  Then  -  =  WN  =  —(*n  +  H'n-i,  and  Wo  =  0.  Also,  note  that  for  Rayleigh  distri- 
a  a  cr 

t=o 

bution,  the  density  function  is  Ja{ol)  =  -y^e_Q  ^  where  =  E(A2).  Changing  the  order  of  integration  as 
illustrated  in  figure  18,  (47)  can  be  rewritten  as 


p=4=  r- r fA(ao)-fA{*N-2){ r 

V27 t  Jo  Jo  Uw*-* 


<t(*-wn_2)  2 

__ ,2/2  /‘  W-i  2qjv_i 

e  /2  /  — ^ — exp( — pt—  ) 


n 


n 


N  — 1 

da/v-i  dx  jdao  dc*i  ♦  •  *  da^-2 

i  /*°°  Z*00 

=  “7 =  /  *••/  /^(ao)*’*/i4(QN-2) 

v27T  2o  2o 

- - v - ' 

TV-1 

drto  doi  •  •  *  dajv-2 


s-*a/»  I  1  - 


exp 


‘  Wn  —  2 


/  <r(z— Wat-2)  \2 

'  *7N-1  ' 


a 


i  z*00  z*00  r°°  2 

=  -=  •••  /  fA{ao)  -fA{(*N-i)\  e~x /2dxda0dai---daN-2 

y/2n  Jo  Jo  Jw’x-i 


N- 1 


j  roc  roo  roc  r  ^2 

--7=1  •••/  f a{<xo)  ■  •  •  f a{.&n -i)  /  exp-—- 

v2tt  Jo  Jo  Jw„-2  L  * 


x2  a2 (x  -  Ws-i) 
VriN-i 


2-1 


dx 


N—\ 

da o  dai  •  da/v_2 


CN 

H 

CN 

b 

/  e_^ 
'o  ! 

r^Cn^ZJ 

N-2 


k—O 

JV-2 


i  i 


^0 


r*’--2 


2  2]]  +2(J2  ^ 


I  ~  H  '4*’ 


(48) 


where 
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Applying  the  fact  that 
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it  follows  that  for  k  >  1, 
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8.2.1  Two  Chips 

For  the  case  of  Ar  =  2,  i.e.,  there  are  two  fading  chips  for  t  £  [0,T],  the  probability  of  error  is 
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From  calculations  on  the  individual  parts, 
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and 


Moreover,  let 
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So,  if  A  <  0,  the  probability  of  error  is 


dx 


if  A  <  0 


if  A  <  0 


f>  _i  _  i  _J}r o _ Wn _ ! _ ijoVn _ 

2  2  y  nr,o  +  2er2  2(r,g  +  rj'j )  s/Shfi  +  2cr2  2(^  4-  vDV^Vo  +  2ct2 

+ _ ®^Vo _ + _ ©V _ 

2A(r,o  +  >7? )  \/2(02  —  A)  2A(ry2  +  r,2) ^2(©2  -  A) ' 

Similar  techniques  can  be  applied  to  cases  with  higher  number  of  fading  chips.  However,  when  there  are 
more  than  two  chips,  the  calculation  of  {A*}  will  be  tedious  and  difficult  to  trace. 


8.3  Analysis  Approach  II 

Another  approach  to  solve  the  problem  in  finding  the  probability  of  error  by  using  the  same  system  model 
described  above  is  to  work  on  the  characteristic  functions  of  the  distributions.  For  the  density  function 
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/^(a)  =  a2/n,  the  characteristic  function  is 

VaUv)  =1  Fi  (l>|?  i^2n)  +  i'02VQe~!L*1 ' 
where  iFi(l,  —a)  is  the  confluent  hyper-geometric  function  satisfying 
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with  r  =  where  rt  =  q^,  the  characteristic  function  of  the  distribution  of  r  can  be  found  by  first 

finding  the  characteristic  functions  for  the  distribution  of  r,.  The  density  function  for  r*  is 


Vi  Vi 


and  the  characteristic  function  is 

$R,(jv)=iF1  ^1,  I;  +  j^-vrifle 

Assuming  independence  among  the  fading  chips,  the  characteristic  function  of  the  distribution  of  f  is 
equal  to  the  product  of  {'P/i,  (jv)},  be., 
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So,  for  r  =  r  4*  n,  it  follows  that 
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which  is  similar  to  the  results  obtained  before.  In  addition,  the  probability  of  error  given  f  is 
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Averaging  over  all  possible  r, 
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Kquation  (49)  can  only  be  evaluated  numerically.  The  procedure  is:  first,  obtain  the  discrete  samples 
of  Vn(jv):  second,  apply  the  DFT  to  get  the  discrete  samples  of  //j(f);  third,  obtain  an  approximation  of 
//j(r)  through  the  Interpolation  []  function  of  Mathematica.  Because  of  this,  it  is  expected  to  have  higher 
calculation  error  at  two  ends  of  interpolation. 
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8.4  Experiments  and  Discussion 

To  compare  the  findings  with  previously  known  results  for  the  one-chip  case,  the  distributions  described  in 
Proakis  [30]  were  used.  However,  since  Proakis  did  not  indicate  the  parameters  of  the  Rayleigh  distribution, 
i.e.,  fi,  there  are  some  difficulties  in  reproducing  his  plots.  On  the  other  hand,  (48)  is  the  same  form  indicated 
in  Proakis,  so  it  is  believed  that  the  results  are  the  same  as  in  Proakis.  Approach  I  was  established  as  the 
reference  and  plots  from  Approach  I  were  used  to  check  the  reliability  of  the  results  obtained  from  Approach 
II.  Another  reference  used  was  the  direct  evaluation  of  (47)  for  single  chip  case.  Figure  19  shows  the  results 
from  the  method  for  a  single  fading  chip  with  fi  =  1  mentioned  above.  It  can  be  seen  that  they  are  the  same. 
For  two  chips  with  equal  fi,  i.e.,  fi  =  1,  Approaches  I  and  II  only  agree  up  to  S/N  =  15  dB.  This  is  shown 
in  figure  20,  where  figure  20(a)  is  for  the  result  from  Approach  I  while  figure  20(b)  is  for  the  result  from 
Approach  II.  Figure  20(c)  shows  both  plots  on  top  of  each  other.  The  deviation  of  the  result  of  Approach 
II  from  Approach  I,  as  mentioned  above,  is  due  to  numerical  and  approximation  error.  It  can  be  observed 
that  this  error  is  more  severe  for  high  signal-to-noise  ratio.  At  present,  it  may  be  assumed  the  above  fact  is 
generalized  for  cases  of  higher  number  of  chins. 

Plots  of  probability  of  error  against  signal-to-noise  ratio  is  presented  in  figures  19  and  20.  Since  the 
equation  of  probability  of  error  were  derived  for  single  and  two-chip  cases,  only  the  results  from  Approach 
II  are  used  to  formulate  the  observations  and  draw  conclusions.  So,  all  plots  presented  in  the  following  are 
obtained  from  Approach  II. 

First,  plots  of  density  function  for  cases  of  different  number  of  chips  is  shown  in  figure  21.  In  figure  21(a), 
{ r;, }  is  set  to  {1}.  In  figure  21(b),  {r/J  is  set  to  {0.5, 0.5}.  In  figure  21(c),  {77*}  is  set  to  {0.25,0.25,0.25,0.25}. 
In  figure  21(d),  {77^}  is  set  to  {0.2, 0.2, 0.2,  0.2, 0.2}.  It  can  be  observed  that  the  distribution  become  more 
and  more  Gaussian  with  a  decreasing  variance.  In  other  words,  with  the  model,  a  very  Fast  Fading  Channel 
is  equivalent  to  an  additive  Gaussian  Channel.  All  these  plots  are  obtained  with  fi  =  1.  Please  note  that 
the  value  of  fi  affects  the  mean  value  r  of  the  distributions.  For  fast  fading,  f  =  1  means  simple  Gaussian 
Channel,  r  <  1  means  a  Gaussian  Channel  with  depreciation  due  to  fading  effect.  However,  f  >  1  means  an 
unlikely  Gaussian  Channel  with  performance  gain  from  fast  fading! 


9  Aspects  on  Optimal  Filters  in  Multiplicative  Fading  Channels 
Without  the  Assumption  of  Synchronization 

The  performance  results  from  wavelet  or  multiscale  signaling  through  a  channel  have  not  been  too  exciting. 
The  results  may  be  summarized  by  observing  that  the  performance  of  multiscale  signaling  is  comparable  to 
t hat  of  other  more  conventional  methods.  From  this  section  onward,  a  more  generalized  problem  contem¬ 
plated  Rather  than  postulating  a  waveform  (such  as  a  wavelet)  and  determining  if  it  performs  well  in  a 
fading  channel,  the  question  of  finding  what  the  optimal  transmitter  and  receiver  signals  are  explored.  In 
the  discussion  here,  some  of  the  basic  investigations  are  presented.  Building  on  the  approach  here,  other 
constraints  can  be  added,  such  as  the  desirability  of  having  the  signal  be  localized  spectrally. 

For  the  channel  model 


r(t)  =  a(t)g(t)  +n(t), 


where  n(t)  is  a  AWGN,  and  where  some  statistic  of  a(f)  is  known,  one  would  like  to  consider  the  problem 

5 

of  designing  a  filter  matched  to  a(t)g(t)y  then  shaping  g(t)  to  obtain  the  maximum  —  ratio  possible.  That 
is,  it  is  desired  to  find 


max 


max 

5(0 

f\g(t)\*dt  =  l 


{max 
MO 

/|MOIa<«=i 


|/Q(*)g(t)/i(0<ft|2 

No 


(50) 


As  “  is  maximized  among  ail  possible  <;(f),  this  optimal  filter  problem  will  also  be  treated  as  a  signal  design 
problem. 
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The  difficulty  involved  for  solving  the  above  optimization  problem  is  that  solutions  must  be  obtained 
to  integral  equations,  which  depend  on  the  channel  statistics.  Instead  of  working  on  the  integral  equations 
directly,  one  can  use  discrete  approximations  to  reduce  the  complexity  of  (50).  As  the  sampling  interval  is 
decreased,  an  asymptotic  solution  is  approached.  The  discrete-time  equivalent  of  (9)  is 

r[/c]  =  a[k]g[k ]  4-  n[k].  (51) 

This  discrete  filter  optimization  problem  is  actually  an  an  eigenvalue  maximization  problem,  which  is 
classified  as  eigenfilter  optimization  problem.  Let  f[k]  =  a[fc]y[A:]  and  let  Rf  be  the  autocorrelation  matrix 
of  /,  and  let  h,  normalized  to  hHh  =  1,  be  the  desired  filter.  Then, 

y[k]  =  h"f  [*]. 

The  above  system  model  is  presented  in  figure  22. 

For  this  system, 

S  hHRfh 
N  "  N0hHh 
hH  Rfh 
N0  ’ 


The  above  Ravleigh  quotient  is  maximized  bv  setting  h  =  xi  with  xi  being  the  eigenvector  of  Rf  corre- 

5 

sponding  to  the  largest  eigenvalue  Ai  such  that  R/X\  =  AiXi.  For  a  given  Rf ,  the  maximum  —  possible  is 

_A]_ 

A  o 

The  optimization  problem  (50)  can  be  approached  by  examining  the  random  process  a[k\  with  autocor¬ 
relation  Ra,  and  finding  the  signaling  pulse  g[k]  with  autocorrelation  Rg  such  that  the  largest  eigenvalue  of 
Rj  -  Rn  o  Rg  is  maximized,  where  A  o  B  =  { ayfty }  is  the  Hadamard  product  of  matrixes  A  =  {ay}  and 
R  =  }.  Kquation(50)  may  be  re-written  to  reflect  this  observation,  i.e., 


max 


max 

9[k] 

£l*Mla  =  i 


f  eig(fla  o  Rg) 

1  No 


(52) 


In  the  above  expression,  eig(.4)  returns  the  largest  eigenvalue  of  A. 

By  obsi^mng  that  autocorrelation  matrixes  are  Toeplitz  matrixes,  and  Hadamard  product  of  Toeplitz 
matrixes  is  also  Toeplitz,  Szego’s  Theorem  can  be  used  to  help  us  solve  the  above  optimization  problem. 

Szego's  Therein  states  that  for  Rn  being  an  mth-order  Toeplitz  matrix,  and  S(uj)  being  the  Fourier 
transform  of  the  coefficients  of  Rn,  and  the  eigenvalues  of  Rn  are  Ani*,  then, 


lim  — 

n-+oc  n 


k^o  Jo 


(53) 


where  /  can  be  any  positive  integer  [74]. 

For  example,  if  r  =  (r„m,  r„m4.i ,  •  •  •  ,rm},  wrhere  =  0  for  JA;|  >  m,  then  the  matrix  Rn  is 
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of  mth-order  and 


"  r0 

ri 

1  o 

r- 2 

r- 1 

T- 2 

r-m 

r —m 

Tm 

^771—  1 

r  m—2 

ro 

r- 1 

. . . 

r  —rn 

Rn  = 

rm 

T  m—1 

rm 

CN  — t 

I  i 

£  £ 

Tm  —  2 

ro 

r-  i 

ro 

*  *  *  r —m 

r_  i  •  •  •  r-m 

(54) 


rm  •••  r\  r0  r_i 

rm  •  •  •  ri  r0  J 

The  Fourier  transform  of  the  coefficients  of  Rn  is 

m 

Sr„  M  =  Y  rke~lku. 

k—  —  m 


(55) 


From  Szego's  Theorem,  the  distribution  of  the  eigenvalues  of  Rn  must  be  very  close  to  that  of  SRn(uj)  as 
the  equality  holds  for  any  positive  integer  l.  In  other  words,  if  entries  of  Rn  may  be  altered,  the  values  of  rk 
should  be  changed  in  a  way  that  maximizes  the  maximum  value  of  S(uj).  Observe  from  (55)  that 


m 

|Sft»|<|  Y  rke-*“\ 

k=~m  (56) 

<  t  W- 

k  —  m 


Now.  consider  iiow  to  find  Rg  for  some  given  Ra  such  that  the  largest  eigenvalue  of  their  Hadamard 
product  is  maximized.  As  both  Rg  and  Ra  are  autocorrelation  matrices,  they  both  are  in  the  form  of  (54). 
It  can  also  be  assumed  that  the  main  diagonal  of  both  matrixes  are  ones,  and  the  off-diagonal  entries  are  with 
values  less  than  one.  Expressing  the  matrices  in  polar  form,  i.e.,  Ra  =  and  Rg  = 

oim*  obtains 

m 

Sro0r,(u)  =  E  (\ok\^\9k\^)e-^, 

k  -  ~m 


which  is  maximized  when  0k  =  —@k  and  \gk\  =  F  At  u  =  0, 

m 

Sr„or,(0)=  Y  lQ‘l- 
k  —  m 

As  all  the  entries  of  Rg  tire  less  or  equal  to  one,  the  above  result  is  the  best  that  can  be  achieved. 

The  following  observations  may  be  made: 

1.  Given  a  fading  channel  for  which  the  magnitude  is  varying  in  time  as  function  a{t)  with  autocorrelation 
matrix  Ra ,  the  signaling  function  g(t)  should  have  an  autocorrelation  matrix  Rg,  which  is  in  the 
conjugate  direction  of  Ra. 

2.  |g*|  =  1  should  be  chosen  to  maximize  the  largest  eigenvalue  of  the  Hadamard  product  Ra°  Rg-  This 
implies  g(t)  has  a  autocorrelation  function  with  constant  magnitude,  which  is  not  very  practical.  On 
the  other  hand,  (56)  indicates  that  there  is  no  benefit  to  require  g(t)  has  an  autocorrelation  with  a 
support  longer  than  that  of  a(t). 


37 


3.  As  some  of  the  <7*  must  be  such  that  \g^\  <  1,  they  should  be  allocated  to  those  small  |a*|s  to  keep 

m 

y;  \ctk9k\  as  large  as  possible. 

k~—m 

4.  As  the  above  guidelines  only  concern  about  the  autocorrelation  function  of  g(t),  which  is  related  to 
its  power  spectrum,  there  is  freedom  to  design  the  signal.  The  above  understanding  of  signal  design 
criteria  will  be  helpful  for  orthogonal  signal  set  design. 

5.  The  optimal  filter  is  the  eigenfilter  corresponding  to  the  largest  eigenvalue  of  matrix  Ra  o  Rg. 


Figure  16:  System  model. 


Figure  17:  Independent  fading  chips. 
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Figure  18:  Change  the  order  of  integration. 


10  Optimal  Signal  Design  for  Multiplicative  Fading  Channels  With 
the  Assumption  of  Perfect  Synchronization 

For  a  multiplicative  channel,  a(f)  (frequency  non-selective)  with  AWGN,  the  received  waveform  when  g(t) 
is  spilt  can  he  presented  as 

r(t)  =  a(t)g(t )  +  n(t),  (57) 

where  n(t)  has  a  two-sided  PSD  of  and  the  autocorrelation  function  of  a(t)  is  known.  If  the  received 
waveform  is  passed  through  a  correlation  detector,  the  output  of  the  detector  consists  of  two  separated  parts 
at  the  instant  of  detection,  i.e.,  the  signal  part  and  the  noise  part,  which  can  be  presented  by  the  following 
equations, 

r8  =  j  a(t)g(t)h(t)dt  and  rn  =  J  n(t)h(t)dt.  (58) 

In  this  section,  some  preliminary  results  are  presented  that  lead  to  the  design  of  optimal  signaling 
waveforms  and  optimal  receiver  waveforms,  when  the  correlation  function  of  a  is  known.  Since  these  results 
are  preliminary,  the  actual  waveforms  resulting  from  these  designs  are  not  presented  (they  axe  the  topic  of 
ongoing  study). 

One  would  like  to  study  how  to  characterize  g(t)  and  h(t)  that  for  given  i?a(r,  A)  =  E(a(r)a(X))j  the 
autocorrelation  function  of  a(f),  the  signal  to  noise  ratio  at  the  receiver  is  maximized.  One  further  requires 
that  g{t)  and  h{t)  to  be  normalized,  i.e., 

G(s)  =  J  s2{t)dt  =  1  and  H{h)  =  J  h2[t)dt  =  1.  (59) 

Maximizing  the  sign al-to- noise  ratio  is  equivalent  to  maximizing  the  ratio  between  the  variance  of  r8  and  rn . 
Because  of  (59), 

eK)  =  y-  (60) 
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Comparing  difference  approach  for  one  chip  case 


(a)  Density  function  for  single  chip  case  (b)  Density  function  for  two-chip  case 


(c)  Density  function  for  four-chip  case  (d)  Density  function  for  five-chip  case 

Figure  21:  Density  functions  for  different  number  of  chips  with  Q  =  1 


Figure  22:  Noisy  signal  to  be  filtered  using  an  eigenfilter  h  [74]. 
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So,  what  one  really  need  to  do  is  to  maximize  the  variance  of  r„  which  can  be  represented  a s  the  following 
functional 


E{r\)  =E{ I J  a(T)g(r)h{T)a(X)g(X)h(\)dTd\ ] 

=  J  j  E[a(T)a(X)}g(T)h(T)g(X)h(X)drdX 

=  j  j  Ra(r,X)g(T)h(T)g(X)h(X)dTdX 
=F(s,h). 

We  use  the  method  of  Lagrange  multipliers  to  maximize  F(s,h),  i.e.,  we  want  to  solve  the  equation 


(61) 


VF(g0,ho)  =  -yi^Gigo)  +  72VH(ho), 


(62) 


where  go,h0  €  L~  (Lebesgue  integrable)  tire  che  extreme  points  in  the  domain  of  F . 
Now,  for  e  £  EL  V  €  C2, 


d  .  .  ~  Fjgo  +  tVJk>)  - 

-Q^F(go,ho)  V  —  diF(go,ho)  ■  V  —1™  f 


(63) 


Furthonnore, 


F(go  +  eV,  ho)  —  F(go ,  ho) 

=  j J  Ra(T,X)(g0  +  eV)(T)ho(r){go  +  eV)(X)ho(X)drdX  -  F(g0,ho) 
=  J f  R«(T,X){tV(T)g0(X)  +  tV(X)g0(T)}ho(T)ho(X)dTdX 
=(  J J[Ra(r,X)  +  Ra(X,T)]go(X)ho(X)h0(T)V(T)dTdX 

=f  J  |y  2/ZQ(r,  A)g0(A)/io(A)dA|  ho(r)V(r)dr 


(64) 


L(r) 


=><3i  F(g0,ho){T)  =  L(r)/io(r). 


Similarb 


=  d2F(go,ho){T)  =  i(r)50(r). 
on 


(65) 


So,  (62)  implies 


L(t)/io(t) 

27i3o(t) 

Mr)go(r) 

27>/io(r) 

|L(r)|  =  2v/|7i72|  and  lso(r)|  = 


|/i0(t)|. 


(66) 

(67) 


(68) 


Because  of  (59),  it  must  be  the  case  that 

|7i|  =  172 1  and  \g0(r)\  -  |/i0(r)|. 

That  means  g{t)  and  h(t)  must  be  similar  to  each  other  up  to  the  same  magnitude  at  the  extrema  of  F(s,  h). 
The  following  observations  may  be  made: 

1. 

L(Xq)  =  0  =>  go(Xo)  —  ho(Xo)  =  0.  (69) 
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2. 


F(go,ho)  =1  J  L(\)g0(\)ho(\)d\ 

=  j  VMglWdx 

=V/|7i72|- 


3.  From  Tricomi  [75,  Chapter  3],  if  RQ  is  symmetric,  then  {<£*},  the  collection  of  eigenfunctions  of  Ra, 
i.e., 


tk*k(r)  =  j  Ra(T,X)d>(\)dX,  (71) 

where  £*  is  the  corresponding  eigenvalue  (distinctive  for  different  k)  which  satisfies  the  orthogonality 
condition 


{4>h,4>k)  =  J  <Ph(r)<t>k(T)dT  =  0  if  h^k.  (72) 

In  addition,  span({()>;fc})  =  £2(a,  6)  with  (a,  b)  being  the  support  of  Ra. 

Now,  let  /(A)  =  go{X)h0{\)  =  £afc0*(A),  then 

L(t)=  j  2Rq(t,  X)f(\)dX 

=  J  2Ra(T,X)Y,a^Wd\  (73) 

=5Za*^^*(r)- 

Because  <■>*  are  orthogonal  and  &  are  distinct,  one  should  be  able  to  characterize  a*  so  that  \L(r)\  = 

constant. 

When  the  above  principle  of  (68)  is  applied  to  a  discrete  equivalent  of  a  multiplicative  fading  channel, 
one  can  find  that  optimal  signals  with  their  corresponding  filters  are  delta  functions.  This  makes  sense 
because  if  there  is  a  variation  with  the  channel,  one  would  like  to  allocate  the  transmission  power  to  a  very 
short  instance  of  time  so  as  to  reduce  the  channel  uncertainty.  This  also  agrees  with  the  finding  of  many 
researchers  [31].  However,  delta  functions  are  not  very  practical  for  communication  engineers.  In  this  case, 
one  need  to  modify  the  cost  function  in  the  Lagrangian  (62)  to  require  the  frequency  span  of  the  signal  g(t) 
being  concentrated  in  a  small  range. 


11  Conclusion  and  Suggestion  for  Further  Study 

In  this  part  of  the  report,  two  methods  of  calculating  error  rate  for  a  multiplicative  fading  channel  by 
assuming  independent  fading  chips  were  presented.  The  channel  statistic  was  also  studied  as  the  signal 
becomes  long  compared  to  the  fading  rate  of  the  channel.  It  was  observed  that  Rayleigh  fading  becomes 
more  and  more  like  a  Gaussian  fading  with  a  decrease  in  variance.  In  other  words,  the  fading  effect  is  more 
and  more  similar  to  an  increase  in  Gaussian  noise. 

By  assuming  the  fading  magnitude  being  of  the  same  sign,  it  was  shown  that  the  optima]  filter  for  any 
fading  channels  are  matched  filters  that  matches  to  the  signaling  function  used.  Moreover,  the  performance 
is  independent  to  the  signaling  function.  However,  when  the  above  assumption  is  lifted,  the  change  of  sign  of 
the  fading  magnitude  will  definitely  make  matched  filters  suboptimal.  As  mentioned  in  [76],  many  researchers 
use  matched  filters  on  fast  fading  problems  believing  that  matched  filters  are  the  optimal  solution  for  this 
situation.  The  derivation  in  this  report  should  be  helpful  to  clarify  the  above  misconception. 

The  correlation  filter  idea  may  be  generalized  to  eigenfilter,  which  is  the  optimum  filter  for  any  random 
signal  when  the  autocorrelation  matrix  of  that  random  signal  is  known  [74].  After  employing  the  notion 
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of  discrete  observables  suggested  in  [76],  it  is  found  that  the  filter  design  problem  is  unified  with  signaling 
function  design  design  problem,  and  they  are  totally  related  to  the  fading  channel  statistic.  This  should  give 
part  of  the  explanation  to  Hansson’s  observation  about  the  probability  on  the  fast  Rayleigh  fading  channel  is 
highly  dependent  on  the  shapes  of  the  modulator  waveforms  [76].  The  design  procedure  for  fading  magnitude 
function  a(t)  with  autocorrelation  matrix  Ra  is: 

1.  for  the  case  of  without  the  synchronization  assumption, 

(a)  determine  g{t)  such  that,  its  autocorrelation  matrix  Rg  is  pointing  at  the  conjugate  direction  of 

m 

Ra ,  with  entries  that  will  keep  ^  as  large  as  possible. 

k=—m 

(b)  design  filter  h  with  value  equal  to  the  eigenvector  of  matrix  Ra  o  Rg  corresponding  to  its  largest 
eigenvalue. 

2.  for  the  case  of  perfect  synchronization, 

(a)  require  g(t)  to  be  localize  in  frequency  domain. 

(b)  construct  a  discrete  equivalent  of  (62)  and  solve  for  optimal  g(t)  and  h(t). 

As  the  above  optimization  criteria  only  involve  the  autocorrelation  function  used  for  a  given  fading 
channel,  there  is  room  to  develop  good  orthogonal  set  of  signaling  waveforms  for  fast  fading  channels.  Of 
course,  studies  should  be  carried  out  about  the  affect  among  the  orthogonal  signal  detectors  because  all  of 
the  signals  will  be  designed  to  match  the  channel  characteristic.  In  other  words,  the  maximum  difference  in 
direction  that  the  different  eigenvectors  point  to  corresponds  to  different  signaling  waveforms  should  be  a 
major  topic  of  study.  This  difference  in  angle  is  critical  in  determine  the  error  rate. 


Appendix  A  Computation  of  the  Variance  Components 

The  fading  component  of  the  matched  filter  output  can  be  written  as 

r(k~q„)T 

r>AV  =  £  -  kT)S(t)pAt  -  iT, -  tD)  ■ 

„=1  i  JkT  (74) 

expL?27r(/Dt  -  ( fc  +  /d)^d)]^,  k  £  28  Z. 

For  purposes  of  computation,  it  is  necessary  to  limit  the  range  of  i  in  the  summation  in  (74)  to  those  symbols 
in  ,s r(t)  which  overlap  with  the  matched  filter  integration  interval  [kT)  (k  +  Q8)T].  If  the  matched  filter  for 
the  fcth  symbol  on  scale  s  aligns  with  the  start  of  a  symbol  on  scale  a,  the  matched  filter  is  said  to  be 
scale -aligned.  Otherwise,  the  matched  filter  is  non-scale- aligned.  A  matched  filter  interval  for  scale  1  which 
is  scale-aligned  with  scale  2  is  shown  in  figure  2(b)  ,  and  a  non-scale-aligned  matched  filter  interval  is  shown 
in  figure  2(c).  If  s  >  cr,  then  the  matched  filter  on  scale  s  is  always  scale-aligned  with  scale  a.  Otherwise, 
the  alignment  depends  on  the  value  of  k:  If  fc  =  0  mod  V* ,  then  the  matched  filter  is  scale- aligned.  We 
also  distinguish  between  those  symbols  that  began  before  the  beginning  of  the  matched  filter  interval  — 
which  are  termed  prior  symbols  —  and  those  that  begin  after  or  at  the  same  time  —  which  are  termed  post 
symbols.  The  timits  of  summation  for  i  in  (74)  are  obtained  by  determination  of  how  many  signals  in  the 
summation  are  prior  to  the  matched  filter  interval  starting  at  fcT  and  how  many  are  post  to  the  matched 
filter  interval. 

Consider  first  the  number  of  prior  and  post  symbols  overlapping  the  matched  filter  when  it  is  scale- 
aligned.  For  a  matched  filter  of  length  q8T ,  there  are  q  —  1  prior  symbols  on  scale  a  intersecting  the  same 
time  as  the  scale  s  signal.  The  number  of  post  symbols  u  on  scale  a  that  overlap  the  matched  filter  is  such 
that 


Qs  >  («  -  1)2", 
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where  u-  1  represents  the  number  of  post  symbols  excluding  the  symbol  with  the  same  starting  time  as  the 
matched  filter.  Since  all  quantities  are  integers,  this  is  equivalent  to 


qs  ~  1  >  (u  -  1)2°. 

The  total  number  of  post  symbols,  including  the  one  with  the  same  starting  time,  is  therefore 


u  = 


Qs  ~  1 
2° 


+  1 


q2*  +  2°  - 


1 


(75) 


where  [x\  is  the  greatest  integer  not  greater  than  x. 

When  the  matched  filter  is  non-scale-aligned,  there  are  q  prior  overlapping  symbols.  The  number  of  post 
symbols  is  determined  as  follows.  Let  i  6  28Z  be  a  number  between  0  and  2°  representing  the  starting  index 
of  the  matched  filter  on  scale  a  relative  to  the  starting  time  of  the  nearest  prior  symbol.  The  number  of 
symbols  u  on  scale  a  that  overlap  the  symbol  of  length  q8  starting  at  i  is  such  that 

q8  +  i  >  u2° 


or 


so  that 


q$  +  i  —  l>  u2a 

q28  +  i  -  1 


u  — 


2 * 


In  the  case  that  i  is  not  in  the  range  from  0  to  2°,  this  may  be  written  as 

q28+(i  mod  2*)  —  1 


u  = 


2a 


(76) 


In  (74)  the  index  of  summation  must  be  changed  so  that  the  index  i  =  0  corresponds  to  the  first  post 
signal.  To  this  end,  we  introduce  the  notation  (k)„  to  represent  the  smallest  integer  >  k  that  is  a  multiple 
of  2'.  For  example, 

(*)o  =  k 


(*)i 


(kh  =  < 


r  k 

if  k  is  even 

{  k  + 

1  if  A:  is  odd 

k 

if  k  =  0  mod  4 

k  +  3 

if  k  =  1  mod  4 

Ar  +  2 

if  k  =  2  mod  4 

k+  1 

if  A:  =  3  mod  4 

k)a‘  and  write  it  as 

^'*  +  1  r(k^)T 

r«.i  (k)  =  Y,  5ZV(*)„.+i  V  /  Ps(t  -  kT)at)Pa(t  -  {k)„'T  -to-  iT,’) 

<7=1  i  JkT 

exp[j27r(/D<  -  (/c  +  JD)tD)\dt 

‘V— 1  rq.T 

=  Y  /  P»(0£(f  +  kT)pa(t  +  (k-  (k)„')T  -  iT„.  -  tD) 

<7  =  1  i 

exp[j2n{fD{t  +  kT )  -  (fc  +  fD)tD)]dt. 


(77) 
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Observe  that 


k  -  (*), 


-J  ° 

-  \  -(2* 


if  k  =  0  mod  2 ° 
(k  mod  2<r))  k  £  0  mod  2a 


so  that  the  effect  of  the  transformation  is  to  shift  the  pa  signal  in  (77)  so  that  when  i  =  0  the  p„  signal  starts 
at  the  same  time  or  after  the  signal  pg;  so  pa  is  a  post  signal  for  i  >  0.  The  index  of  summation  in  (77)  can 
now  be  changed  to  reflect  the  number  of  prior  and  post  signals: 


rim 


N.+l  .g.T 

(k)  =  ^2  &«,,(*)„,+«*'  I  P*{t)Z(t  +  kTfy'rit  +  (k  -  (k),,')!  -  iTa>  -tD) 

<r=l  i=lk(tr')  0 

exp[j27r(/£>(t  +  kT)  -  ( fc  +  /o)t£>)]dt 


(78) 


where 


J  -(q  -  1)  k  =  0  mod  2C  (scale-aligned) 
k'-a’  ~  |  _g  Jt  ^  0  mod  2”  (non-scale-aligned) 


and,  from  (75)  and  (76), 


«s,*(<7) 


f  |2^=lJ 


/c  =  0  mod  2^  (scale-aligned) 


2,^^(<rmod2',)  _i  j  ^  q  mod  2CT  (non-scale-aligned) 


It  is  straightforward  to  show  that  r,$(k)  =  r,ti(k  —  2®  ). 

The  autocorrelations  are  computed  conditioned  upon  a  known  set  of  interfering  bits,  which  is  the 
set  of  all  bits  referenced  in  (78).  The  correlation  between  fading  components  can  now  be  written  as 

r,(k:  1, 1|XM)  =  ^E[rl1(h)rt'1(k)\lt,k] 

N,  +  l  AT#  +  1  u»' 

=  £  £  E  E  6 

<Tl  =  1  <72=  1  ii—lk{(7\)  i7=lh(<7,<2) 

[  [  p'a{ti)p»{t-i)Ri{<2  ~  ti)p’ai(ti  +  {k  -  (k)^)T  -  iiT^  -  tp) 

Jo  Jo 

PoAt-2  +  (k-  {k)^)T  -  i2T -  tD) exp(2j7r /o(t'2  -  ti)]dt,cit2. 

The  correlation  between  fading  components  at  lag  one  can  be  written  as 

r ,(k.  1,2| I..*)  =  rl(k; 2, 1|Z,,*)  =  \e[t]a  (k)r,.2(k)\I,,k} 


(79) 


S.*-l  ft'.+  IV.tki)-1  U.M»-3*'|(<T2) 

=  E  E  E  E 

ai~l  <72  =  1  |l=/fc(<7;) 


*Vi ,(*)„<  +ii2‘T>  ft<r2,(*-2*')„i+i22'’i 


m.t  fv-T p.{ti)Mti)Ri(t2  _t]  _2 >'T)p'c,x{U+{k-{k)<)T -i,T^-tD) 

Jo  Jo 


(80) 


P.,(t2  +  ((*  -2 »')-(*-  2®'),- )T  -  t2T,.  -  tD) 
exp[2j7r/o(t2  -  fi  -  2®  T)]dfi<ff2 

Appendix  B  Programs  Used  to  Evaluate  the  Probability  of  Error 
for  Fast  Rician  Channel  of  Section  3 

The  probability  of  error,  presented  by  equation  (19),  when  applying  multiscale  wavelet  signaling  is  evaluated 
by  R1 1 .12 .  c  and  comp.pe.Bit .  c  which  appear  in  algorithm  1 .  Due  to  the  tremendous  computation  involved, 
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many  tricks  are  used  to  speed  up  the  computation.  First,  Rll_12.c  is  used  to  calculate  (79)  and  (80) 
without  specifying  Z**.  The  result  from  Rll_12.  c  is  stored  in  a  table  corresponding  to  slot  of  different  shifts 
of  overlapping  that  will  be  imported  to  comp^>e_Bit .  c .  The  program  comp_pe_Bit .  c  will  take  the  table 
output  by  Rll_12.c,  assigns  different  set  of  Z8  *  and  completes  the  calculation.  As  the  size  of  ZS)*  can  be 
huge,  it  is  impractical  to  calculate  the  probability  of  error  for  all  the  different  combinations  in  Z*,*.  So, 
what  comp_pe_Bit .  c  really  does  is  randomly  generate  a  large  enough  number  of  elements  of  X8 and  use  the 
average  probability  of  error  obtained  from  this  subset  of  Zaj*  to  approximate  the  true  probability  of  error. 


Algorithm  1  Compute  probability  of  error  for  Rician  fading 


/< 


•  Program :  Btl_12.c  —  sod  If  lad  fro*  Dr  Nooa’i  *aioa2.c 

•  Chat  Lo 

•  Utah  Stata  Oalrarslty 

•  Data:  Oct  14,  1*96 

•  Tbli  program  eoapnta*  tlia  autocorralatioa*  upon  a  toon  tat 

•  of  iatarfarlag  bits  and  output  to  a  tabla  that  till  bu  mad 

•  by  coap.pa.81t . c  which  alii  amigo  to  aach  • lot  of  tha  tabla 

•  a  ra'.atlva  data  aalaa  (-1  or  1)  to  cacslata  tha  probability 

•  of  arror  by  raadoaly  gaoarata  a  larga  no* bar  of  trlala. 

. / 


•lac.ada  <(a»t.tb  l> 
flat. ad*  <a.1bh> 
•lac.ada  <d imdtypaa . h> 
•lac.ada  <atdlo.h> 
•lac.ada  <afd’.lb.h> 
■lac.ada  <aath.h> 
•lac.ada  citrtng  h> 


•daft aa  ca..oc. array  pea Lioe. array 
Maflaa  cm. .oc.aatrll  peal  1  oc .matrix 

Maflaa  aiala.b)  ((a)  <  (b)  T  (a)  :  (b)) 
Maflaa  aatla.b)  ((a)  <  (b)  *  (b)  :  (a) ) 
•dafta*  tPS  1  Da-10 
Mafia*  U'IC.I  10OO 

/•  |.oba.  aarlab.aa  •  / 
tat  0. 

tat  li 

tat  I. pat. pat . 

char  aataaaaa  [  )IG  J . 

(at  aaacoaf , 

doab  a  pal  .pkl  .  taaap.  ad.t 
doab.a  'paldata. 
doab.a  >paidata. 
char  nattaap  t  *C)  . 
char  oatf  l  ,a[.V)J  . 

F  lf  -~p*p. 

F  ’  t.f  »oatfp 


tat  R.tltab  .ail a*  a  lOOOO, 


doab.a  •R.tltab. a.  /•  fadlag  corral  at loa  valaaa  •/ 

/••For  t.tl  tab. a  oa.y  ••/ 
doab . a  t.ti.aca.a. 

. . . . / 


doab . a  baa a )0 ■ doab .  a ;  , 
void  ialt.A.tltab.aldoab.aJ 
doab.a  ft. tl (doab.a) . 
lat  .  oaar  1  lat  ,  til), 
lat  appardat,  lat.  lat!. 
lat  h.ataai.balgllat ,lat J . 
void  lait.pat.pbl (char  •  ). 
/•••  aaaa; at  faactloaa  •••/ 
doab: a  pal (doab.a. , 
doab.a  pal (doab.a) , 
doab.a  pat  I (doab.a) ; 
doubla  phi  1 (doab.a) ; 
doab ! a  pa!3(doab.a) ; 
doabi*  phi  2 (doab.a) ; 
doab : a  •caia.pat; 
doab) a  acala.pall ; 
doab) a  ac*.a_p»12, 
doab.a  acala.phl; 
doab: a  acala.phl l ; 
doab.a  acala.phl 2; 

/ . . 


lot  pr laa ( lat )  ;  /•  aagaa  '■(■agM<IaTsag»h:Ii)  •/ 

doabla  T.iagaadat);  /•  T_<»ag«ML>»2‘ {»ag»a}T  •/ 

doabla  d#: tat ; 

doabla  qgaaa(doabla  (*faac) (doabla) ,  doabla,  doabla); 


/•  la  It  lata  t.tl  tabla  •/ 

/•  rat  an  tala*  for  R.rt ,  fading  corraiattoa  •/ 
/•  \ ( 1 _h ( aagaa) \)  •/ 

/•  \(a.{*  ,k) ( aagaa) \)  •/ 

/•  \( (•- (b) .Vtagaa) \)  •/ 


/•  Tbay  ara  atactly  tba  iiat,  •/ 

/•  dip  11  cat ad  to  tpaad  ap  iatagratloa  •/ 


/•  iangth  of  aaTalat  algaal  •/ 

/*  latarayabol  tltaa  for  basic  scala,  ahaa  TO-I, 

latartyabol  tltaa  for  scala  1  sill  ba  2  •/ 

/•  no* bar  of  «cala>.  It  •/ 

/•  so* bar  of  polsts  for  psldata/phldata  •/ 

/•  atm  bar  of  coaff  for  ptl/phl  •/ 

/•  aaaplad  tiaa  for  pal/phl  alth  T«t  •/ 

/•  aavalat  foactloa  •/ 

/•  acaliag  faactioa  •/ 
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TOid  gauleg(*Old) ; 

double  r_poi  [GiULSGJtf  .  aeg[GKJL8Gjr3  ; 

. . . . * . 

The**  global  rariablea  are  for  the  eultl  Integration 


double  il,i2,Bl,B2,Cl,C2.Tl,T2,T3; 
lot  in_* ,  in.aegeal,  la_aeg»a2; 
doable  tart; 

double  (•  P_t)  (double); 
doable  (•  P.tageal)  (double); 
double  (•  P_aegea2)  (doable); 
double  fiat  111 (double) ; 
doable  f let 112 (double) ; 
double  flntl21 (double) ; 
doable  flat  122 (double) ; 
let  BlgT; 

double  rntoep(lut,iut,lut,iattiBt,int);/»  rlicoaqj(kta,aageaI , aegea2.il .12)  •/ 
double  rl2coap(lBt,iBt.lBt,iBttlBt,iBt) ;/•  rt2eoep(k,a .aegeal , tegeel. 11.12)  •/ 

. . . 


iat  eela(lat  urge,  char  •  •argr) { 

char  •ptlfaaee; 

tat  ptlaaeechar .type; 

lat  ■  .k.aegeal  .  aegaa2.it  .12; 

doable  Bd ; 

let  couBt"0: 

char  • 1 l#Be»e[20j ; 

lat  id jle. tat ; 

tat  t; 


Iflargc  —  1)  { 

priatf (“uaage:  111. 12. c  peraatf llenaeeVa”) ; 
et it  tO) ; 


ptifueee  -  (char  •)e«Teed(arga[l]  ,~palfBaae“,tp8luaaechar  .CHlB)  ; 

••reed  arg* l I ] , “70“ , tTO , DOUBLE ) ; 

*>r*wl .  argi  [  I J  ,  "ll"  .Ul .  (BT) ; 

e treed  ergr  [  1 )  ,  ”Bd“  ,tBd  ,  XTJBU )  ; 

at reed  erg*[l] . "type"  .rtype, (IT) ; 

aareed  erge [ l ] . .ta , (IT) ; 

a treed : argr( l] . "De.teT” .hdeltet .OOUBLZ) ; 

aareed  '.argali  J  .  ~id J3el  tet“  .Bid JOeltat  , IIT)  ; 

/Malt la. lae  eoaclaaat  aid  aelghtt  of  the  Geaai-Legeudre  quadrature  foreule*/ 
gee. eg 

tf , hdjle.tet : 

de. tut  »e  prtee(t)/?  0), 

•tge~.tegee  priee  ai); 


latt.pal.pui  pt l feeee 1 
telt  .a  titeb  e  hd  . 

tpriatf  *t  .eeae»."l»hdtg-.ar(r»tlj  ,»d)  ; 
tt  r  c  py  .  oe*  Me.*!.  Meet .  . 

•  treat  oet»l  e.~  oat*:. 

•trept  oeteep  ,*t  .«e*  . 

•treat  oe teep . *  eep": 

eap* p  •  * ■'pM  <elatp,'t'l  . 

•hi  e  ’eep'p  • 

eap*p  •  *opea  oeteep ,  . 

pr  l tt  *  ~'~*>+  d  eot  opee  eapf  1  .a\t“  5  ; 

*pr  i  at  *  eep*fl» .  'O-tg.  at»ld .  t«Xd\a“  ,  eereaaea ,  TO  ,1*  ,t)  ; 
fc./vee  eep»p 

oet*p  •  *opee  oetM.e.-e". 

oettp  •  topee  oat  *  1 . e . . 
prletf  “Coe.d  eot  opee  oetf t . e\e“) ; 

ft . ote : oet  *p  . 


If'type-t- 

for:i*-l  i<  itt'i*  prieetJ J ) .!*♦)< 
r-t*.  M<priee  « 
eeptp  •  *opee  oeteep. ’a";, 

•ht.e  ’eepfp  < 
eepfp  •  topee  oeteep. "a") . 
prletf f "Coe . d  eot  opee  eepfl.eVa”); 

r 

fprtatf  eepfp. -*-ld\f.i)  ; 
fc.ote (eepfp. . 

for ( tegee 1 • 1 .  tegeel  <■(■••»  l ) , tegeel **){ 
eepfp  *  fopee  oeteep .  "a  M , 
ehl.e! 'eepfp;i 

eepfp  *  f opee. oeteep. "e"  i 

priatf ("Coe. d  eot  opee  eepf 1 ;e\a") ; 

> 

f  priatf  (eepfp,  Magee  !«Xd\t"  .tegeel) ; 
fclote(aepfp) ; 

for ( aegeaT* 1 : aegee2<«(Ba«l ) ; tegea2*+){ 
eepfp  ■  f opee. oeteep. “e" ) ; 
uhllaC 'eepfp! \ 

eepfp  »  fopee (oeteep. "e") : 
priatf ("Coeld  hot  opea  eepf l ie\a") ; 

> 

f priatf (eepfp . “tegea?»ld\t " .segea?) ; 
fc(oae(eepfp) ; 


?  or  (11*lo«*r  (prlM  (angMl) ,k) ;il<uppnr(pri»n(i#gMl) ,prlM(a) ,k);ll+*){ 

«4pfp  ■  fopnn(ootMp.'V')  ; 

■kiln  (•  Mpfp)  { 

**P?P  -  f  opnn  (oatBBp ,  "»")  ; 
print? ("Could  not  opnn  MpfilnW); 

> 

fprlnt?  (Mpfp,  "il*Xd\n'*,il) ; 
felon#  (Mpfp) ; 

for(12*lo»nr  (priM  (»#gM2)  ,k) ;12<npp*r(pri»n(*ng**2)  ,prl»n(i) ,k) ;12**){ 

Mpfp  »  f  opnn  Coutmnp , "»") ; 

«hlln(  imapfpK 
Mpfp  «  ? opnn ( oo tm»p , "n")  ; 
print? ("Could  not  opnn  »*pfiln\n”); 

> 

f  print? (Mpf  p,  "12*Id  :",12) ; 
felon# (Mpfp) ; 

print? ("XgVn"  ,r llco«p(k ,• , mgMl  ,  1 1 ,12) ) ; 

print? ("count*Xd  \n ", ♦♦count ) ; 

> 

*tp?p  *  f opnn (outMp,"**')  ; 

»blln(  ?Mpfp){ 

Mpfp  *  fopno  (oatmnp  ; 

print? (“Could  not  op*n  MpfilnW); 

> 

fprlnt?  (Mpfp,"\n") ; 
f  clonn(Mpfp)  ; 
ontfp  ■  fopno  (ootf  11*  ,•'**•) ; 
nbll*£  'oatfpH 

ootfp  ■  fopnnContfll  ; 

print? ("Conld  not  opnn  ont?iln\n“) ; 

> 

fprlnt? (ontfp, "\n“) ; 
felon# (outfp) ; 

> 

> 

> 

#.  »*' 

for!l*Q,1<(|<< (In-prlM (n) ) ) ; 1**) { 
k»1*<  l<<prlM(«) ) ; 

Mpfp  ■  f op*n (ostaap , "n") ; 

(Mp?p)  ( 

Mpfp  ■  t opnn 1 on t map ,  "a"  1 ; 
prlntffCon.d  not  op«D  MpfilnW); 

f  pr  l*t?  '.Mpfp ,  ”k*Xd\t” ,k) ; 

? t . on* I >*p?p) . 

for  :  ■•(m  1*1  ;  nngul  <*(ln«1)  ;nagMl«*){ 

Mpfp  ■  f^pMlNtnip,***); 
vfei <«*p?p; t 

Mpfp  •  fop**  lontmap  ,•*•*’)  ; 

print? ; “Co*. d  not  op*n  MpfllnVi”); 

fprlnt?  !Mp*p.  “tnjfM  1  *Xd\t“  ,  nngnui  l )  ; 
ft . onn ;awp?p . . 

for  '  ukm>  1  .  i«pn]«!inM)  ;n#gM2*+){ 

Mp'p  ■  *»p*n .ontntp ,“n“) ; 
nil  'Mpfp  ■; 

Mpfp  ■  fop**  :  ontMp  1  ; 

pr  1  nt ? : “Con . 4  not  op**  MpftlnXa**); 

fprlnt?:  Mpfp  .  ~nngM2*Xd\t  “  ,  nngnnlJ  ; 

? t ,  on#  Mp'p .  . 

for  '  it-  onnrprlM  n*gMi )  ,k)  ;  1 1  <*ppnr  (prlnt(i«gMl )  ,prlM  (a )  ,k)  ;  il**){ 

Mpfp  •  f  op*n  ;  »n tmap  .'■*”>  , 
nm  'Mp'p  «. 

Mpfp  •  fopnn  o#tMp.“n") : 

print? '“Con. a  not  op*n  Mpfl.n\n“); 

f  print?  Mpfp  .“t  l-ldW.ll)  . 
ft.o**iMpfp : . 

for  1 1  2* .  o**r  :  prtM  :  ••gnnl!  .  (k-  ( 1  <<pr  1m  ( ») ) ) ) ;  t2<npp«r (priM ( mgnn3)  ,  print*  (■)  ,  £k- ( 1<<prtM  ( i ) ) ) )  ;  12*+H 

Mpfp  •  fop**  MitMp,  “*“) , 

nni.n: >*tp'p:< 

Mpfp  ■  fop**  ontMf  .“n'3  . 
print?!  “Con  .  d  not  op**  Mp?i.n\n“J. 

fprlnt?  Mpfp.  “12-Id  “.121; 

fc. on*: Mpfp  •  . 

print*:  "XjW.rl  Jtoop'fe.n.nngMl  .  n*gM2 ,11,12)1  ; 
print?  “connt -Id  Vn'.Mconnt!, 

> 

Mpfp  •  f  op*n  :<y*tMp.“n“l  ; 
nhi  .n;  'Mpfp :  ■ 

Mpfp  •  fop**  .  o*tMp,  “n") ; 
print*! “Con. 4  not  op**  Mp?l.#\*“); 

y 

fprlnt?  ;Mpfp,“\*“) , 

f  c .  onn  :  Mpfp :  . 

ontfp  •  fop*a (ontfi ; 

■*ll*(  'ontfp}.; 

ontfp  ■  fop** loot fl , 

pr  lit  f  l"Con .  d  not  opan  oot*  1 .  n\n" )  ; 

fprlnt? (ontfp, “\t“3 , 

(ontfp} ; 

> 

> 

> 

> 

> 

Mpfp  -  fop*n  (ontMp .  )  ; 

•btlaj  !Mpfp){ 

Mpfp  ■  f opnn (catmap ,  "*")  ; 
printf ("Could  not  opnn  Mpfll*\n”l; 

> 
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f prlatf  (mapf p .  "DOTH ! !  • !  !\u'*)  : 
fcloaa(Hopfp)  ! 


void  lnit.R.xltabla (doobla  BdK 

/.  initiate  tha  B.xi  tibia.  Tha  tabla  b».  R.xitablaalaa,  and  lmltlatad 
According  to  tha  iu  tao  raqnlrad  (  \(q_«  T  \)  ahara  a  la  tha  hat 
bcaU).  that  la,  R.xi[0j  for  R_xl(0).  B.xi U.xltablaal**~l)«R_xl (-aitao)  . 

•/ 

1st  1; 

doobla  tao.total.tiJkA; 

total.tlaa  *  (doobla)  (TO*(l«Ia)*Q)+BlgT; 
tan  ■  total.ti»a/(B_xltablAal*A~l) ; 

R.xitabla  -  (doobla  •)c*lloc_arTay(MroBLK.R_iltablaal*a,,'a_xltabla“) ; 
ford  ■  0;  1  <  R.xitablaaioa;  !-*-►)  { 

B.xltabla [1]  •  baaa jO(2«R.PI*Bd»i*tao) ; 

> 

B.xi.acala  *  (R.xitablaalxa-l ) /total. ti»a ; 


doobla  B.xi  (doobla  ta«H 

/.  iaataa  tha  foactioa  la  liaaar  bataaan  aaighborlng  data 

B.xi (tao)  »  B_xltabla[V  ♦  dalta* (B_xitabla[nvl]-R.ritabla[nJ ) 

•bar a  a  •  lat  part  of  tao*acala,  dalta  »  fraction  part  of  tan«acala. 


lat  n ; 

doobla  dalta.  BCalad_tl»a; 


tao  afaba(tai) ; 

acalad.tiaa  »  (tao*R_xl_acala) ; 

■  •  (lat!  acalad.tlma; 
if  lo  —  'B.iltablailaa-l) ) 
rat ora (R_*l t abla [oj ) : 
a. »a  if ;a<{B.*ltablaalaa-1)){ 
data  •  aca.ad.ttHa  -  a; 

rat  on  CR.xltablalaJ  *  dalta* (R_xltablaC«+t]-R_litrbla[a] )) ; 


prior? '.”11. tl  oot  of  ranga!o»Xd  tao«Xg\oM  ,a  .tan)  ; 
rotam  0; 


lat  . oaar : tot  tag**,  lat  k) 
lat  as. 

■  a  a  :i«iago»; ; 

oh 1 .  a  tk  <  0 .  k  ••  o3 ; 
k  -  kill. 

i?;"kl  ratoro  -50-1'; 
a.aa  ratoro  -0. 


lat  oppar 1 1  at  tag** ,  lot  a.  lat  k) 


tot  »»|o». 

U«gH  " 

ohi.ak  <  0  toagoa . 

k  •  klttagaa 

i  * :  •  k . 

ratoro  . lot : : : ;0<<o:  •  t tag**- 1 5 /t oagma) ; 


lot  k.alaoa.katg: lot  k.tot  tagnal 
/•  cnotpota  k  -  Ik;. tag*#  •/ 

l 

lat  o3. 


■  3  •  i”«*go*. 

oht.aCk  <  o;  k  **  t3. 
k  •  k  t  03. 
lfCk)  ratoro  0. 
ratoro tk  -  a3J . 


void  lolt.pol.phi5 char 


H 


r:LI  •la.pkl.pbl.ft.a; 
lot  1; 

cbar  otr [tOOj . 


/•  polflla  format: 

■  ■  i Oort  taaa  of  aavalat  flla 

a  a  anabtr  of  coafflclaato  (lotagar) 

pi  i  pOl.ioap.ad.t  •  taMplad  tljoa  (doobla) 

0  ■  langth  of  aaoalat  (la  lotagar*)  (lotagar) 
l.pal.phl  ■  oombar  of  oasplao  (lotagar) 

■  tap  la  I ,  oa*pla2.  ...  (doabla)  (oaralat) 
a asp la 1 .  »aspla2.  ...  (doobla)  (ocAling) 

•/ 

/•  tha  data  la  aada  aalag  tha  program  “sakaaava”  •  / 


ia.pal.phl.flla  ■  fopaa(f llaaa»a.“r") ; 


lf(in_pai_phi_fila  —  IULL)  { 

printf ("Could  not  op«D  ln.pai.phi.f ila\n”) ; 

•tit (0) ; 

> 

ficuf  (ln.pfl.phl.f  11«>“X*\b>>  ,ur«iM)  : 

fgeta(atr,100,ln_pal_pbi_flle) ;  nuacoef  »  atoi(atr); 

fgeta («tr , 100 , in.pal.phi.f 11*) ;  pal.phl.aeapled.t  ■  atof(atr); 

fgata (atr , 100,lB_pal_phl_fll«) ;  0  ■  atol(atr); 

fgata (atr, 100, lu.pal_phi.flla) ;  I.pal.phl  «  atoi(atr); 

paldata  •  (doabla  •Jealloc.array (DOUBLE, (*_pai_phi*l) ."paid at a") ; 

pbldata  ■  (doubla  *)ealloc_array( DOUBLE, (I.pal.phl*!) ,"phidateM) ; 

for(l  »  0;  1  <  I.pai.pbi;  i*-0  { 

If ( ! fgata (atr , 1 00, in.pal.phi.f ila) )  ( 
prlntf ("Cannot  raad  »a*pla\n") ; 

•xlt (0) ; 

> 

paldataCl]  *  atof(atr); 

> 

paldata tl.pai.pbl]  «  0; 

ford  »  0;  1  <  I_pal_phl ;  l-*-*)  { 

if (!fgata (atr, 100, ln.pal.phi.fila))  { 
printf ("Cannot  raad  aaaplaW) ; 

•tit (0)  ; 

> 

pbldata[lj  ■  atof(atr); 

> 

pbldatatl.pai.phl]  *  0; 
fcloaa(ln_pal_phl_f  ila) ; 


doabla  pal (doable  t){ 

/•  laaowa  tba  function  la  llnaar  bataaan  Balgbboriag  data 

pal (tea)  ■  paldata[n]  *  delta* (p*idete[n*l]~p«idata[n]) 

akara  a  *  lat  part  of  tan*acala,  dalta  ■  fraction  part  of  ta«*acala. 

•/ 

lat  n; 

doob.a  dalta,  tcalad.t 1m ; 

•  ca.ad.ttM  ■  t/acala.pal; 
a  •  (let)  acalad.tUaa; 

If (■>■!. pel. pbl) 

rat  om  1  p»  1  data  [E  _pa  1  _pb  1 J  )  ; 

a.aa  * 

da. ta  ■  acalad.tUaa  -  a; 

ratam  (paldataUJ  ♦  delta* (paldata (a* I J -paldata [n] )) ; 

> 


doab.a  pki: doabla  tl{ 

/•  liiiM  tba  function  la  1 laaar  bataaan  neighboring  data 

pbl (ten)  ■  pbldata[a]  •  dalta* (phldata(n* 1] -phldatatn] ) 

•  bare  a  ■  int  part  of  tao*acala,  dalta  «  fraction  part  of  taiMcala. 

•  / 

lat  a. 

doab.a  dalta,  acalad.tlM; 

acalad.tlM  *  t/atala.pbl; 
a  ■  Hat!  iu.td.tiM, 

If (•»■«_ pal .phi ! 

rat  am  (pbldetall.pil  .pbl  j )  ; 

a.aa  ( 

da . ta  ■  aca.ad.ttM  -  ■; 

rat ara'phldatataj  ♦  da.ta* ( pb tdata la* ij -pbldata [aj ) ) ; 


doab.a  p«11 (doab.a  tit 
lat  a. 

doab.a  da.ta.  tea . ad. t tea . 

Kt.ad.tlM  *  t/tca.a.ptl  I  . 
a  ■  datJ  tta.ad.tiaa, 
if (a>*R.p«i .pat : 
ret«m(pa1detell.pal.pbl  j  ;  . 

a.aa  t 

dalta  •  tca.ad.tlaa  -  a. 

rat  am  ( p(  Idat  a  La  J  •  da . la* ( ptldatela* 1 J -p« Idata [aj) ) J 


doab.a  pail (doable  t ) < 
lat  a; 

doab.a  dalta.  aca.ad.tlM; 

■calad.tlaa  »  t/acala.ph 1 t ; 
a  ■  (latl  acalad.tlM, 

If (a>al.pai.pbl) 

rat  art (pb 1datall.pil.pbl J J , 

a.aa  ( 

dalta  ■  tct.ed.tlM  -  a; 

ratam  (pbldata  (.a)  *  delta*  (phidata  (a*  I J -pbldata  [a] )  )  ; 

> 

> 

doabla  pal2(donbla  t)( 
lat  a; 

doabla  dalta.  acalad.tlM; 


acalad.tlM  ■  t/ecaia.pal3; 
a  ■  (let)  acalad.tlM; 
if (•>■■. pal.pbl) 


retnrn(paidate£l_pnl_plilJ)  J 

el»e{ 

delta  -  »caled_tl»n  -  *5 

rtton (pnldututn]  ♦  deltn*(paldutu[:i>-M]-p»iduta[ii])) ; 

> 


doable  phi2 (double  t){ 
lat  •; 

doable  delta,  •ceJ.«d_tl»e; 

■celed.tlM  ■  t/ucnla_pbl2; 
a  ■  (let)  tceled.tlJM; 

If (B>»I„p«l_pbl) 

return (pbldete [K.pa l.pbl] ) ; 

else{ 

delta  ■  ecnled.tlJM  -  n; 

return (pnidete [a]  *  delta*  (phidutu[o-*l]-pbldutu[i»J)) ; 

> 

> 

lat  prise (1st  negseH 

/•  negsa ’ •( »eg»n<J«T»eg»e: aegse- 1 )  •/ 

lftiegse  <  It) 
retnrn  negsn; 
elae 

return  I* ; 


double  T.iegneltut  negseH 
/•  T.(neg»e>«2<negs*}T  •/ 

return  (TO* ( t<<neg»n) ) ; 


donb.e  f tat  I22(dooble  teuH 
doeb.e  Integral; 

letegre.  ■  *_*1 (teu-BtgT) ; 
tutegra.  *eP.i ( t *art*tau) ) ; 

Integra.  •■F.»epea2 ( v»arteteu*T2) ) ; 
return  Integral; 


doeb.e  flutist ; double  *H 
doea.e  Integral ; 

•  art  •  t . 

Integra.  ■  qgae* (flat  ITS, Cll-t) , (B2-t) ) ; 
tutegra.  **F _• ! t  J ; 

Integra.  *»P.negse1 ( (t*Tl J 1 ; 
return  Integra. , 


doeb.e  ? let  t  iS'.doeb.e  caul1, 
doeb . e  I • t egra . . 

Integra-  ■  N.tlltaul, 

Integra  **r  . •  :  eert *taul  . 

Integra.  **r . eegsaS I  I  tart *tan*T2) ) ; 
retnrn  Integra. 


doeb.e  ft a* ill. doeb.e  t .  *. 
d»nb  e  letegrn . 

nett  *  t 

Integra  •  qgeet  1  flat  I  1 2 , 1 II  -t  1  ,  )  )  ; 

tutegra  •«*.!  t  . 

let egra.  >•*  eegsel : : t • ' 1 ) ;  . 

return  Integra. . 


doeb.e  rlkAUg!|*t  n.lnt  n.  lat  negset.  Int  HgaaS .  Int  11,  lut  12){ 
/•  r  I \ctmf  ,  a  .  •  .  i«fH  i  .  negjaeJ .  1 1  ,  IS)  •/ 
doeb.e  Integra. . 
doe be  nte.elmg. factor . 

If  (»<■**: I 
F.n  •  pel. 
le.n  ■  n. 

> 

el  net 

P.n  *  pat . 

> 

ucele.pnl  ■  (donb.e)  (1<<te.»). 
tce.e.pnl  •■pnl.pbt.nesp.ed.t ; 
tea. a. pat  •  non. a. pit. 

If leegee! <■■•)< 

P.negse l  ■  pat! ; 


> 

el*e{ 

P.negM i  ■  pan  . 

la.negBai  *  la; 

> 

ecale.pnll  ■  (doable)  ( I«la_negnw1 )  ; 
aca.le.pnll  *«pn  l.pbl  .ie«p  .ed.t ; 
tcoae.pall  •  ■CAle.pnll; 


If (eegna2<«B«){ 

P.ligMl  “  pal2; 

la_aegne2  *  aegna 2; 

} 

elae{ 

P.Mgail  «  phl2; 
ln_iegne2  *  Ba; 

> 

acale_pei2  ■  (doable)  (l<<in_aagne2) ; 

■cala_pal2  •■*pai_phl..«enpled_t; 
acale_pbl2  *  ecale_pal2; 

T1  ■  TO»k_niaaB..kBlg(k,iB..eegnet)-it*T_eagna(in_eegnel) ; 

T2  ■  TO*k_niBBa_kalg(k , in_aegne2) •*12»T_eegne ( la.iagatS) ; 

At  ■  nax((-Tl).0); 

12  *  ■la((TO»(tJ«lB_aegMl)-TI),(T0»(O«iB_a)))i 
6 1  «  nax((-T2) ,0) ; 

B2  -  Bln  ( (TO* (Q<<ln_»egna2)  -T2)  ,  (T0«  (Q<<in_e)  )  )  ; 
integral  >  qgem  (f  lnt  11 1 ,  11.12); 

/•  For  a* a<3  with  alnpleO*/ 

/• 

iategral  ••daltat ; 

Integral  ••daltat; 

•/ 

/• 

Tkla  takaa  care  tbe  ac-iliag  of  pal  end  pbl  >blch  waa  not  done 
la  pat  0  and  phlO  . 

•  / 

a  c* .  el  ag_  factor  e  (doable)  (  l«(lo.fegnal«la_eegna2)  )  ; 
ace. a lag. factor  ■  aqrt (acelelng.f actor) ; 
acalelag.f actor  •«  (double) (1<<la. a) ; 

Integral  /*  acelalag_f actor ; 

napfp  ■  fopao (catnap, "a") ; 
aht.a( 'napf p) { 

napfp  ■  fop«e(oatmap,lae") ; 
prietf (“Coeld  not  open  napf  1  ie\n’")  ; 

fpriatf (napfp, “Xg  ".letagral); 

oatfp  ■  f opaatoatf  lle,“a‘*) ; 

•hi .a; 'o«tfpi< 

oat f p  •  fopaa (oet? tie, “a") ; 

pr t at f  !  "Coo . d  eot  open  ootf lle\o”) ; 

fpriatf ; oat fp,"Xg  " . letagral ) ; 
f t . one (oatf p ; 
rat  an  Integral  ; 


doab.a  ruenapdat  k.tat  a.  lat  aagnal,  let  aagna2,  lnt  11,  let  12){ 
/•  r12eonp(». a. aagnal  ,»agna2,l  l  ,12)  •/ 

doab . a  tatagra. . 
doab.a  tea. a  1 ag. factor ; 


t  ■  •  pa i  , 


ace. a.  pa  l  »  . doab.a; 

tea. a  p*  I  *«pt 1 _ pa t _ tanp . ed_ t , 

ate. a. pat  * 

If; aegna 1<-Ba: • 

P_ eagna  <  •pall, 
la. aagnal  •  aagnal 


P. aagnal  •  pail, 
la. aagnal  ■  . 

> 

ace. e. pall  -  (doeb.e!  (1«le. aagnal)  ; 
aca.a.palt  •■pa  1 .pa  1 . tanp , ad. t , 
tea. a. pail  ■  ace. a. pall. 

tf (aegna 

P.aegna3  •  p»12, 
la_aegna7  *  aagnal, 


P.*egne2  •  pa  111, 
la.aegn*3  •  la. 

■ceie_pat2  •  (doab;a)  ( l<<  la.tegnn2) ; 
ecela.pal2  •■■?  at.phl.eanpled.t ; 

•ca.a.pbl2  ■  acela.pt 12; 

T1  ■  T0*k.n1aaa.ka1g(k, ta.aegnet l-ll-T.aegnada.aegnal) ; 

T2  •  T3«k.alaai.k«tg(  (k-U«te.*))  ,le.aagne3) -12*T.aagna(la.»agne2) ; 

11  ■  nax ( ( -T 1 ) ,0) ; 

12  •  nl a( (TO* (0<< ln.iagna  1 > - T 1 ) ,(T0» (Q<< ta_«) ) ) ; 

•1  ■  nax ( ( -T2) ,0) ; 

12  •  nln((T0»(0<<la.*agne2)-T2) , (TO- ((K<la_a) ) ) ; 
letagral  •  qgaaa (flatlSl ,  11,12); 

•cal el og_f actor  ■  (doable) ( 1<< ( la.aagna I *la_aegne2) ) ; 
a cel el ag. fee tor  •  kqrt (tee) elog.f actor ) ; 
teal al Dg_ factor  •■  (doable) ( 1<<1b.b) ; 
letagral  /*  ecelelag_f actor ; 


WXpfp  ■  Jop«B(0Bt*«p.,*a")  ; 

■hilat !*apfp){ 

»ap?p  ■  f  opaa  (ootmap ,  "a'*)  ; 
priotf {"Could  not  opaB  »*pf  HaNn”) ; 

> 

fpriatf (»apfp,"Xg  ".Integral) ; 

f  c.loaa(mapfp) ; 

ootfp  ■  fopafi(ootfila,"*") : 

»bila(!oDtfp){ 

ootfp  ■  fop*B(oatf ll«p“4") : 
priatf ("Could  Bot  opai  oBtfilaW) ; 

> 

(printf (oBtfp."Xg  ".iutagral) ; 
fclosa(ootfp) ; 
ratara  intagral; 


•  Inatrlctl  intagratlon 


told  gaalag(vold) 

< 


lat  J.i; 

doabla  *1.  s,  xi,  pp.  p3,  p2,  pi; 

« ;  ■  l ; 

for(i«0;  t<CiIUC_*;i*-*){ 

i  ■  coa(>».F>(l*0.'’,6)  A2»GiOLBG_I-K>.&)) ; 

4o( 

pl«l  0. 

p2~0  0. 

forlJ-O  .  j<«2-GiCL*C_*);J**H 

P^P3. 

p-'-pl 

p1«:.2  0*j*l>a«p2-)-pJ)/(J-*l>; 

pp*3»GiL'LXG_«»Cx«pl-p2)/tx»*-l)  ; 

•  -  ** . 

x-x.-pl/pp. 

•d I . a> fab*  tx-xl ) >KPS) ; 

t .  pot  [Gi-1  *a  .1-1-  I  J  ■  T>«; 

Bag  JC.i-1- lj<2. 0«il/{(1 . O-X-x) -pp»pp) ; 


tioit . ■  qgaaa;doab.a  S •fame)  (doabla) ,  doabla  a,  doobla  b) 


lat  J. 
doab .  a 


xa.  a*,  a. 


wO  *•  b*a:  . 
rr-C  *•  b-a.  . 

•  •0  . 

for  ’  JH5.  JfClJlJC.l 
a*-Tr..  po.t Jj 

a  *»  a*gi}]»  '  •faac  .  x»~d¥ )  •  ( *faac)  (x»-d*) ) ; 
ratara  *  **tt 


deab .a  a*aa - •  a; 


*oab  •  tx.f ,aa* .aaal .a»a2. 

if  : . •!«*•»•  i . :  «  l  o'  < 

?•*•» 

0*y«  <  -  113«2i903&4  .  0*y*  (6&1«1»640 . 7 
+?•'-  113144  24  l**y*(T73«  13017*y.(-|!H.»0S24£«)))))  I 

*ai2-fcT*A»«»04l 1  b*y«  ( 103*W2*«>  b*y«(M»4e»0.7t8 

♦y'kvrr?  64«fc>»7«(MT  .tua?!  3*7*1  .0))))  ; 

aaa**aa i/aa«2 . 

•  a  »*  ■ 

**«  0/ax 

r*** 

IT^I-O 

aaat-1  0-y-:-0  HXW*i2a«27»-2.*7.  (0  7734M0407*-* 

*y«:  -0  207jjT(M3»a-b*y0.2093e872na-«))); 

tail  ■  HD  lb*;M»»w»fca-1»7«(0  .  H304S87&&#-3 

»j«-0  l HT661#-fr*7* (0 . 7®2J0*fcld to-« 

•f0  »J4»3ili2a-7)J)j 

a.  •••qrt : 0  &3M IPTT2/ ax) • (cot (rx) •aaa1-*«ala(rx) *aaa2) ; 
ratara  ui . 

> 


/• 

Local  Varlab.aa 

coaplla-t  :wu4  “apcc  -o  1111.12  -g  Kit. 12. e  -I/atr/laclada/local  -lpaatllb  - 
; coapt 1 a-co^aaad : “gcc  -o  HI  1.13  -g  R11.I2.C  -I/afiViaclada/local  -L/opt/llb 
; coaptla-coaaaad  “gcc  -o  dKII.12  -g  R11.I2.C  -I/iBr/laclado/local  -L/opt/lib 
End: 

•/ 


11a  -1*“ 

'IpMtltb  -11*  -1*" 
-Ifmat.lb  -11*  -1*" 
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/' 


•  Program:  comp_pa_Blt . c  —  modlrtad  from  Or  Moo*’*  mat on 2. c 

•  Chat  Lo 

•  Utah  Stata  UnlVaralty 

•  Data:  Oct  14.  1W8 


•tacloda  cpmatllb . b> 
•incloda  <allb.h> 
Vinclnda  <aiadtypai . h> 
■lnclnda  <*tdto.b> 
•incloda  <«tdllb.h> 
•lnclnda  <math.h> 
•lnclnda  <»trlng.b> 


Maflna  calloc.array  pcalloc. array 
Maflna  cal loc_matrlx  pcalloc.matrlx 

Marina  mlB(a.b)  ((a)  <  (b:  ?  (a)  :  (b)) 
Marina  max(a.b)  ((a)  <  (b)  T  (b)  :  (a}) 
Maflna  BPS  1.0a-10 
Maflna  GiUUKJI  1000 


/•  gob*,  tarlabiai  •/ 

tat  0;  /•  langth  of  varalat  algnal  •/ 

doabla  TO:  /•  lntarnymbol  tlma  for  baalc  acala,  than  TO  1 . 

latartymbol  t Ima  for  acala  1  will  ba  2  •/ 
doab.a  10, 
doab .a  id 
doab.a  S. 
doab.a  3. 

tat  /•  aombar  of  acalaa.  Is  •/ 


lat  . oaar t lat ,  tat) ; 

tat  appar.lat,  tat,  tat); 

tat  b.aiaai.balgdnt.tat}; 

lat  kt  tg ; tat ,  tat}; 

lat  pr  lx  :  t  at  J  , 

doab.a  T . aagma  flat!; 


/•  \Cl_K»agma)\)  •/ 
f •  \(..{i,k}(«ga.)\)  •/ 

/•  \( (k-(k)_\aag*a)\)  •/ 

/•  aagma  '  ■  ( aagmadiTaagma  : Is)  •/ 

/•  7_{aagma}«2-(aagma>T  •/ 


doab.a  . .  /•  rll[»]  [lj  [aagmalj  [up«2]  •/ 

doab.a  ••••••r13; 

sold  latt.rlli'fear  •!, 

*otd  latt.rlilchar  •), 


doab.a  » .» ,» .pit  tag. a. at a. 3a. taPht ; 
doab.a  laoab.a.lat .lit ! . 
doab.a  flat  .doab  , 
doab.a  da. tat  ■  I  a- 3. 


doab.a  ym  doab  a  .  *ftit!  (doab.a)  ,  doabla,  doabla}  ; 
void  gaa . ag ' told  . 

doab.a  l.poalCiClK.I]  ,  aag iGiUUtG.I]  ; 


doab.a  Ml.ta.oa  tat.  tat  , 
doab.a  rl2.ta.aa: lat .lat ) , 
doab.a  r22.  ta.  tat ,  tat ; . 
doab.a  comp .Pa ' told : . 
doab.a  craap  .O.  tarsi  '  tat  ,  tat}; 

totd  • at  bitmap ' s«td : .  /•  randomly  aat  tba  bitmap  tabla  •/ 

/••  ait  ttaf*  ••/ 

tat  ••bitmap.  /•  bit  pattarna  •/ 

lat  •mai.poatt  toa  ,  •mu .  potl  t  loa  ,  /•  for  cal  tba  dlmaaatoa  of  bitmap  •/ 

totd  tatt. bitmap  told; . 


told  1 a 1 t . b I tmap  told  < 

tat  a.t.t.MgMi  .tagma2.11.t3. 

1st  tamp l .  tamp 2 . 

(•  Utad  to  ttora  tba  max/mta  poattloa  for  lndltldaal  acalaa  •/ 
tat  j. 


max. poattloa  •  (tat  • Jma. loe (a lamof ( lat } • (laM ) } - 1 ; 
ala.  poattloa  •  (tat  •  ;mal  loc  Catxaof  (1st )  •  (IiM) ) -1 ; 

for (J»l ; j<*(la*t) ; J*a){ 
max  _  poa 1 1 loa ( J  J  •  0. 

■la.posl t loalj J  •  0; 

) 

for  (a-l ;  «<• (I»*l )  ;■«-*)  i 

ford  —  1  ;!<{!«  (la-pr tma (a )  )  ) ;  1  — ) { 
tM»(  1«prtma!a)  3 ; 

for  ( tagma  1  ■  I  ;  a agma  I  <»  (la«  1 )  ;aagma  !«-•){ 
for (aagma2-l ; aagma2<a(la*t} ;  sagma2*“» )  { 

for (t l»loaar( primal tagmal ) ,1) ; 1 icnppar (prlma(aagmal) .prima(a) .1) ;t1*a}{ 
for(12*lobar(prlma(aagma2}  .1) ;  U<nppar(prlma( aagma2)  .prlma(a) ,k) ;  !>•-»){ 
tamp  1 ■ (ktlg (k , pr tma ( tagma! } } *1 1 • ( l<<pri*a( aagmal ) )  )/< l<<prlma( tagma 1 } ) ; 
tamp3-(ka  lg  (k  .prtma  ( aagmaU} }  *12.  ( t«prtma  (aagma2)  ) )/(  l«prtma(  aagmaH) }  ; 
If (tampicmln. poattloa [sagmal]) 
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ain.pot It lon[*«g»&i3-teapl ; 

•111  If  (t«*pt>mM_poiitioB[BBgMl]) 
au.potitlOB  [aagMl]  »t«apl ; 

If  (taatp2<*in_poaltioB  [«agM2]  ) 

■1b  _pot  it  loa  [aagM2]  »t«p2 ; 

alia  If  Ctaatp2>max_p<>BitioB[aagM2]  ) 

■ax.poait  Iob  [lagMa]  -xmmp? ; 

> 

> 

> 

> 

> 

> 


for(a«t ;■<* (!■+!) ;a-*-*H 

for(i»0;  i<(l«(Ia-prl*a(a) )) 
k*i*(l«priM(i)) : 

f»r(Mgul*l :  aagMl<»(ia-H)  ;««pul^H 
for(kagM2-l  ;  «agM2<-(If»l)  ;tagM2+-*H 

for(ll*lo»*r (prla«(*«gMBl) ,k) ;il<«ppar(pri»a(BagMl) ,prl»«(*) ,*) ;ll-»-0{ 

f or C t2« loaar ( pr  1m C i«gM2 ),  (k- ( 1  «pritaa (■))));  12<opp«r ( priaa ( aagM2) , priw (■) .  (*- <  1  «pri»a (■)))>: *2*+H 

t«apl-<kalg(k,prl»a(iagMl))*il*(»«prlM<BagMl)))/(l«priM(aagMl>); 

taaap?*  (kilg ( (k-{  t«priM(a)  ) )  .priM (i«gM3) ) *12* ( l«priM (i«gu2) ) )/ ( 1  «prl»a C»agp»2) ) ; 

If  ( t  «p  1  <*la_po  *  1 1  loa  [  I*gM  1  ]  ) 

■ia_po»ltion[*«g**lj*to«pl ; 

•  Iia  If  CtOBpl>«AX_pelltiOB  [t*gHl]) 

MX  _poa  1 t 1  OB  [  ■  «gM  1 J  *t  ««p  1 ; 

If  (tMp2<«iB_po»ltloB  r*  jgmAl]  ) 

■1 b _po • 1 t i OB l ■•£■■23 "t ««p2 

alaa  If  Ct*atp2>MX_po*ltlo«t»«g»»23  ) 

MX.potltlo*  [*«g»»2]»t*p»; 


> 


prtatf(“Id  Xd\n*'  .■in.poattlon [a]  ,Mi_po«it lot  [i] )  ; 


bitMf  »  Clat  ••)m*  loe(»l*»of  Ctit 
f  or  :  *•  1  ;  a<«CIaM  )  ;  a*-»)  i 

DitMp[aJ-;iaf)M::oc(«i*aof  (lBt)»(MX_poBitioB[a]-miB_poaltioa[t]  +  l))~*iB_po.ltioB[a]; 


1st  »agM,  iBt  k) 

i 

tat  a2 . 


*3  •  t  ; 

wkt.aCk  <  0;  k  *-  «2; 
k  •  H«2. 

tfi'k)  ratara  -(0-1;: 
•  ratara  -0 . 


tat  appar :  tat  i«cm.  tat  a,  lat  k) 

tat  tup*, 
tup*  •  H«up*  . 
akt  .a.  k  <  0.  k*-  t  up*, 
k  ■  tltiagM. 

tf ' 1  k 

ratera  1  tat  ;::!0<<a  :  ♦  ti.ga*-l)/tMp»); 


ratara  :  tat :  :  ( 10<<a  :  •  k  -  D/tiapa); 


tat  k.Btaai_baig:tat  k.tat  aag»a; 
/•  coapata  k  -  Ik'.aagM  •/ 

1 

tat  a2 . 

•  3  •  I  «»apa  . 

Bkt'.aCt  <  o;  k  *•  a). 

k  -  k  X  a3. 

tf(»k)  ratara  0. 
ratara; »  -  *3,  , 


tat  kalgMat  t.  tat  up*;< 
tat  top. 


taap  •  k.ataaa.kftgtk  ,i*p*J  . 
rasp  k; 
ratara  C-ta«p) ; 

> 

lat  prlaallat  aagM.'i 

/•  a*pt  ‘  ■  ( aagM<la7aagM  :  tagM- 1 )  •/ 
ratara  tala  (  up*  .la  ) )  . 

> 

doabla  ?_aagM(tat  i«p*){ 

/•  T.(a«pa*<J'(«*p*)T  •/ 

ratara  (T0»  ( t«iagM) ) ; 

) 

void  tatt.rt  l  ( ckax  •»aj*a){ 
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lot  a, k, aagaal  ,taga*2.il  ,12; 
lot  1; 

eh  or  filanaaa[30]; 

FI LB  *infila; 

rll  ■(doabla******)aalloc(alxaof (doobla*****) •  (■•■*!) )-l; 
for  ( a*  1  ;•<■(!••»  1)  ;a*-*H 

rll[i]  "(doobla*****)aalloc(alJt#of  (doabla****)*  ((l«(la-prl*a(a))  )*l)  )*1 : 

/•  -M  ao  that  th«  raallait  ralaa  for  1  cao  bo  >1  •/ 
for(i*-l;K(l«(Ii-prlaa(a)))  ;l*-*){ 
k*l*(l«prl*a(a)) ; 

rll [a] [i]*(doabla****)«alloe(al*aof (doobla***) •(■■♦l)}-!; 
for (aagaal*l ; aagaal <*(I»+1) ;i«ga4l+*){ 
rl l[a] Cl] [aagaal ]* (doobla* **)aalloe(a laaof (doobla**) *(Ia*l))-t j 
for(aaga*2*l ;  aaga^c-da+t)  ;  «aga*2*-»)  ■{ 

rll [a] [1] [aagaal] [aagaa2]*(doabla**)aalloc{al*aof (doobla*) • (appar (prlaa(aagaai) ,prlaa(a) ,k) -loaar (pri»a (aagaal) ,k) ) )-loaar (prlata(  aagaal) ,k) ; 
for (11 "loaar ( pr loo ( ■ ago* 1 ) ,k) ; 1 1< appar (prtaa (aagaal) ,pri»a(a) ,k) ;il*-*M 

rll [a] [i] [itpul] [aaga*2]  [llJ*(doubla*)aalloc(aL»aof (doobla)*(oppar(priaa(aagaa2) ,prlaa(a) ,k) -loaar (prlaa(aagma2) ,k)))-lowar(prlJBa(»ag*a2) ,k) ; 

> 

> 

> 

> 


for (a*l ;a<»(Ia*l) ;■*-♦){ 

aprlotf  (f  ilanaaa  .''IiXdBdXg"  ,1«M  ,*  ,Bd)  ; 

•treat (f Ilanaaa," .oot”) ; 
lof  1  Io*f  opon  (f Ilanaaa  ,  "r")  : 

If (lofllo— *ULL){ 

prlntf  ( "Sot  abla  to  opoa  X*\n" ,f Ilanaaa) ; 
axit(O) ; 

> 

for  U*- 1  ;l<(l<<(Ia-prtaa(a)))  ;1*-*H 
k*l*(  l«prlaa(a)) ; 

for <aagaa1»l ; aagaa 1 <■ (■•*! ) : aagaa  !*-*){ 
for(aagaa2*l  ;  aagaa2<-(la*l )  ;  iagaa2*-»)  { 

for ( 11 -loaar (prlaa( aagaal) ,k) ; 1 l<oppor (prlaa( aagaal ) .priaaCa) ,k) ;ll*o){ 
for(12*loaar(prl*a(tag»a2)  ,k)  ;  12<nppar(prt»a(aag»a2)  , prlaa(a)  ,k)  ;12+*){ 
facoof  (inf  ila,"XlaM ,krll  [a]  [1]  [aagaal]  [i«p»2]  [ll]  [12]) ; 


> 

f c.oaal lof 1 la) ; 


void  talt_rl2(  cOor  ,ua«){ 

lat  ■  . k  .  itpu  1  ,  aagaa? , 1 1  . 12  ; 

1st  1; 

char  f 1 . aaaaa [30  J ; 

HLB  •loft.o; 

rl2  ■ (doab: a* •••••) mI loc( a laaof (doobla* • ••») •(■•♦!) )-t j 
for  ( ■*!  ;•<*(■•«!)  ;•■**){ 

r!2[a]  “(doab .a»**** )aal  loc  (  a  laaof  (doobla****)* ( 1«  (Ia~prlaa(a ) )  ))  ; 
for(l»0;l<(l<< (Sa-prlaa (a)));l*-»){ 
k-l*( l <  <pr laa (•)) . 

r!3laj ll J*!doob. a** **)aal loc (ataaof ( doobla***) •(*•+!) )-l ; 
for  (*ag»a!*l  , aagaal <*(■** 1 )  ;  aagaal*-*)  { 
rt3laj[tj [aagaa ! J-(doabla***)aal loc  t  * laaof (doab la** )•(!•*  1) )-l ; 
for  !  aagaa?*!  .  ■afp<}<>!li«l )  ,  aagaa?*-*)  { 

rl2(a]  t  tj  l aagaal )  [ aagaa? j  ■  !  doob.  a**)aal  loc  (a  laaof  (doobla*)  •  (appar  (prtaa  (aagaal ) .prtaa (a)  ,k)  -loaar  (prtaa(  aagaal )  ,k)  ) )  -loaar  Cprlaa  (aagaal)  ,k) ; 
for i lia.oaar I prlaaJ aagaal ) ,fe) ; 1 1 < appar (prlaa( aagaal ) .prtaa (a) ,k) 

r  1 2 [ a 3  llj  [aagaal  J  [aagaa?]  [  1 1  )-£do«bla>)aalloc(alaaof  (doobla) •  (appar  (prtaa  (aagaa2)  .prtaa (a)  .  (k- (l<<prlaa( a) ) ) ) -loaar  Cprlaa (aagaa2) ,  (k-( l«prl*a(a)  )) ) ) ) -loaar  (prlaa (aagaal)  ,  (k-  C  l<<prlaa  (a) ) ) ) ; 


> 

> 

for ! »-l ; •<•;■••! 1 ;a»* ) • 

aprlatf : f 1 . aaaai , "laid  Wtg" .taaa.a .  kd ) ; 

•treat (ft. aaaaa . "  oat" 1 ; 
laf 1 ;  **f  opaa { f 1 .aaaaa . *r“) . 

If  {taf  l'.a— rJU  )( 

prtatf("lot  ab.a  to  opaa  Ia\a** , f  1 1  aaaaa)  ; 
at  1 1 ( 0  J  , 

> 

for  ( 1-0;  l<  ( !<< (la -prtaa  !  a) ) )  ;  l*-*)< 
k-l*{ 1<<prtaa(a) 1 , 

for  (aagaa t*1  .  aagaa  1  <•  (la*  t )  ;aagaa1«){ 
f  or  (  Mgu>  1  ;  aagaa?<*  («a*t )  .  aagaa?** ){ 

for ( 11  ■ loaar  (prtaa ( aagaa 1 ) ,k) ; 11 coppart prlaa (aagaal) . pr laa ( a ) ,k) ; 11*-*)  { 

for (12* loaar ( pr laa C aagaa2 ) , (k-(i<<prlaa(a }))) ; 12< appar (prl»a(aagaa2) . prlaa ( a) , (k- ( !<<pr laa ( a) ) ) ) ; 12^-») { 
f  acaaf  (laf  l.a,"X.a".krl3[a]  [1]  [aagaal  J  [aagM2]  [1 1]  [12]  )  ; 

> 

> 

> 

> 

> 

fcloaadaf  1  la)  ; 

> 

> 


doab  la  rl  l_ra_laa(lat  a,  lat  1)< 
lot  aagaal,  aaga^2,  ll,13,k; 
doab  . a  taap,  a «.■ 11 0  ; 


k*i • ( l<<prlaa (a) ) ; 

for  (aagaal*!  ;  aagaa  1  <■  (Ia*l )  ;  aagaal*-*)  { 
for (aaga*2-l;aagaa2<- (»•*!)  ;aagaa2-**){ 
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fordl-lowarCprlmaCaagaal)  ,k)  ;il<nppar(prl*a(Bagm*l)  ,pri*a(a)  ,k)  ;11+*H 
for(i2»lo»ar(prima(fagma2)  ,k) ; 12<nppar(prima(aagma2)  ,prl*a(a)  ,k)  ;12-*-*){ 
t«p  »  rl 1 [a]  [i]  [aagmal]  [i*gu2]  [11]  [12]  : 

taap  .-bitmapdBgmalJUkaiglk.primaCtagmaimiW^prlmaCaagaaim/CK^rlmaCaagmal))]: 

taap  •■bitmap  [i«pu2]  [  (kalg (k.prlma(8agm*2)  ( l<<prlma(aagma2) ) )  /( t«prima  (aagm*2)  )]  ; 

>“  +“  tmmp; 

> 

> 

> 

> 

ana  **D : 
ratum  an; 

> 

doobla  rl2_aaloa(liit  a,  lot  1){ 

1st  aigMl,  *ag**2,  11,12,*: 
doabla  tamp,  anm-O; 

*■1*  (l«prima(a) ) ; 

forCaagmal-l :»agaal<*(Ia+l)  ;aagmal-*-0{ 
for  C««pml2*1  ; tagma2<«(Ia*l) ;  aagm*2+*H 

for (1 t-loaar (prima(aag»al )  ,k)  ;il<oppBr(pri*a(f8gmal)  ,priao(a)  ,k) 
for(  i2-j.oaar(prlma  (aagm*2) ,  (t-( l<<prlma(a)  )) ) :  12<nppar  (prlma (sagma2)  ,prima(a)  ,  (b- ( l«prl*a(a) )) )  ;  12++H 
tap  ■  r 12[»]  [1]  [aagmal]  [aagm*2]  [11]  [12]  ; 

tap  •  ■bitmap  [aagmatj  [CtatgCt  ,prl*a( aagmal)  )*ll*  ( l«prl*a(aagm*n  )  }/(l«prlaa(aagmal)  )  ]  ; 

taap  — blt*ap[iagma2][CkBlgUk-U«priaa(a))).pri»a(Bagm*2))*12«(l«pri*a(Bagma2)))/Cl«priaa(aagma2))]; 

■  a  taap; 

} 

> 

> 

> 

> wa  »«0; 
ratara  am; 

> 

doob.a  r22.*alaa(lat  a.  lat  k){ 
ratara  rl t.aalaaCa, (b-l)) ; 


Boat . a  floats. a  iocalpal,  lat  a,  lat  t) 

Boat. a  btl.k22.kl2; 
doab . a  rbol ,  rho2; 
doab.a  latagral; 

bit  -  rll. am. aa ! a , b  J «I0 ; 
b 1 2  ■  rlS.ralaat a .b J ; 
b22  •  fO-»r22.Ta-laa(a  ,k)  ; 

ir;tn<a; 

bit  •  -bit; 

it ; b22<o : 
a 22  •  -b72. 


rbo t  •  S/ktl. 
rboJ  •  S  Ak  22 . 

ata  •  1 1 2/  aqrt 111  I )*aort (b22) ) ; 

2  •  a*  r*o 1 *rbo2 ; . 

*  ■  *•  ra«2-rbo t . . 

V  •  s^rt ’ rbol »ra«I 1 . 
pataag  a  •  .oca. pat. 

latagra.  •  <ypaa«  flat,  •■.H/3,  *  _P1 /25  / 1 4*b_PI )  ; 
ratara  latagT*. . 


doab.a  Mat  Ooab.a  t 

doab.a  I 

doab .a  tl ,  raaa  t 

/•  tiiatt  that  ia.tarbt  •  0  a*4  pataag.a  •  o-pl/2  •  / 

/•  for  }*.uni  •  »-pt.  pat  aag.a— -•-pl/3 ,  cat U>al taPbl-palaagla}-©*/ 
■  a  -  v«ata - t i 1 . 

If  li—O 

prtatf  '  -b^\a"  :  . 

tt  -  V*alaC>a.taPti-paiaag:a)/l; 
tl  *•  ;«ta*lti  Ipalaag.al ! . 
raaa .  t  ■  arp  -a  *t  1  . 
rat ara; raaa. t) . 


doab.a  coap.^a.taPbt ' tat  a.  tat  k){ 
lat  taap; 

taap  •  bitmap [a] It] -MtmaplaJ  [b-t] ; 
tf  !taap— 01 
ratara  0. 
a.aa 

ratvra  B_P: ; 


doabla  eoap.Patrotd J { 
lat  a.k; 

doabla  Pal.  P*”0 .  taap; 

forlt»l ;  a<-(ia-»l)  ;a«-»)  { 

Pat-0; 

for(b-0;b<(  1 << («a-prlma  !  ») ) )  ;b**H 
DaltaPfcl  ■coap.Daltai’hl (a.k) ; 
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t «p  -  F((W_PI/2),s.k)  -  F((-K_PI/2),«,k); 
Pal**feba(ta*p) ; 

> 

Pel  /-  (double)  (l«(Ia-prl»e(a))) ; 

Pel  /■  (doable)  (l«pri»e(f)) ; 

Pe  *«  Pel; 


retain  Pe; 


void  set bitmap (void H 
lat  i,k; 


for(a*l ;•<■(!•  +  !) ; e-»-0{ 

for(k^io_po*ltloo[«]  ;k<eBez_posltloa [»]  ;k«-*H 
lf(r*ad(X(lUlD.HiI/2)) 
blt«ep[t] [k]  ■  1; 

else 

biteep [■] [k]  ■  -l; 

> 

> 

prlatf ("\n") ; 

for  (>■!  ;§<*(!•-»  1)  ;«♦*){ 

for (k^ln.posltloa [•]  ;k<-«ax_posltlon[s)  ;k-*-*){ 
prlBtf("Xd  ".bltmeptaj LkJ) ; 

> 


prtotf  ("\b**}  ;ff  laah(stdott)  : 
> 


»ele;  1st  erge,  ckar  ••argv){ 
doab.e  SMb; 

doable  Hart  BblO  .endEblO , stepEblO ; 

doab.e  db.Rb.Pe; 

tat  rwaaetber; 

tat  leaberof rno,  eoant»0; 

doab.e  •oatPetable; 

lat  t. oatPetable; 

tsar  •rllfoaae,  •rl2fnaae; 
cser  peraef  t  .eUO)  ; 

9  :  LI  •oatPe, 

1st  1. 


Iflargt  —  1)  { 

prlatf ("asa^e:  coap.pa.Btt  par  oaf lle\a“)  ; 
at  It :o,  . 


•trepy (peraaf 1 .e.argrtl ] ) ; 

rltfsaae  •  user  ■) arread (argil tj .“r 11 f otsa" .krnaaeekar, CEAfl) ; 

ri  2-faaae  •  Ukar  • ) erreed (ergv 1 1 J ,"r 12f itse" .krsaaeehar .CHIP) ; 

as  reed  per  serf  1.  a.  “M"  .  tSd  ,  DOUBLE !  ; 

asreed . paraaf  l.e/'ls"  ,fcla .  IIT) ; 

eareed  paraaf  t  .e,‘‘Q~,kQ.  I  FT); 

e«  reed : paraaf I . a .“TO” . iTO , DOUBLE) ; 

eareed  '  paraaf 1  .  a ,  “S“ ,tS .  DOUBLE)  ; 

a* read  pereaif  1  .a ,  "SOdb"  .ESDdb  .DOUBLE)  ; 

eareed  .  pereeif  t . a  ,  "at  ertfhSO"  .ft start EblO  .  DOUBLE) ; 

eareed  paraaf  1 .a , "aadEbSO" .teed EblO, DOUBLE) ; 

eareed  paraaf t . e . "at epEbiO" , ftatapEblO .DOUBLE) ; 

eareed  peraaf 1 . a . "Inaberof rae~ , tJeabarof raa , I IT) ; 

/•tel  tie.  Isa  abattaaaa  ead  salgbtt  of  tka  Geaaa-Lagaadra  qaadretsra  fomli 

•« 

lalt.rn  ri  ifaasa;  . 
tatt.ru  . rlDfaaaal  . 

Ib  •  S.  /••CHECK  •/ 

3  ■  S*po# : (0 . -SDdb/tO; . 
t  at  t . bt  teept i . 

■  pr  l  at f ' per aaf  1 . a . "la Bdl^SDIglsld" . ergv 1 1 ] . Bd . SOdb . Is ) ; 

■treat (per ear l.e.~  btt  oat"), 

l.oatPeteb.a  •  ( 1st i ( oedlbSO-etartthEO) /atapBbEO  ♦  t; 
oetPeteb.a  *  ( doeb . e*  ,'ea ..»<(•  lseof ! doab a) •  l.oetPetable) ; 
f or ( 1«0 . 1 <■ . oetPat ab .a, 1 *• ) 
oetPatab.altj-O. 
do*. 

set bt reap ( ) . 

for (db  •  tiartBblO.l-O,  db  <■  eedBbJO,  db  *•  atapEbEO , lee)  { 

10  •  Eb»poa (  10 . ( -db/ 10) ) , 

Pe  •  coap.Pe ( ) ; 

oetPetableliJ  •  oatPeteble[lJ •  Pe; 
prt.tfrX*  Xg\a”  ,db  ,Pe)  ; 

> 

coas  t  *• ; 

oatPe  •  fopeatperaaf  He, "a")  ; 

If  (oatPeeerULLH 

prlstf!"Eot  able  to  opea  Xa\n" .paraaf 1 la) ; 
a*tt(0); 

> 

fprtatf  (oatPe,  “\**coest>ld  Bdlg  E»td\M\a" .coaat ,Bd .!•) ; 
for(  teO:  1  Of.  oat  Pet  eb.  a;  !♦■*) 

f prist f ( oatPe ,"Xj  Xg\n",  ( start BbEOM  •atapBbEO) . (oot Pat abla(lj /coast) ) ; 
f close (oatPe) ; 

)ak  I  1  a  (  coant  daaberof  mi )  ; 


> 


h 


•  Btt*arieal  iatagration  (froai  Xraarical  Racipaa) 


«/ 


Told  ganlag(TOld) 

{ 


Ut  j,i; 

doobla  *1,  x,  xl,  pp,  p3 ,  p2,  pi; 
xl  ■  1; 

1  or  ( 1-0 ;  KGiOLBG.I ;  1~) { 

X  -  eo« ( H.PI •(!+<). 76) /(2*CiOLK_l+C. 6)); 

6oi 

pl»l .0; 
p2-0.0; 

for(J*0; J<(2*G10LEG_*) ;  J*-*H 

p3-p2; 

p2-pl; 

pl*( (2.0»J*l)*x»p2-J»p3)/(J*l) ; 

> 

pp*2*GAGL£G_B» (x*pl-p2) /(a^x-l) ; 

xl-x : 

x»xi-pl/pp; 

}»blla(faba(x-xl)>BPS) ; 
i_poalGiSLEl_f-l-l]»  xi'x; 

aaglGi:.’L*G.I-l-l]»2.0»xl/((1.0-*«*)*pp-pp) ; 


doab.a  qgaaa  Cdoabla  (•fxBc) (doablx) ,  doabla  a.  doablx  b) 

i 

u»  ). 

doab.a  tt,  SB.  dx.  •; 


*■>•0  b»;b*x), 

rr»0  *• : b-a; , 

•-0. 

for  J-O.  J  <G  I  .l-ES  .1 ,)«-»)  { 

dr»TT«X . pot t j]  ; 

■  *•  aag[  j  J  •  :  :»fxmc!  (x»»dx) ♦ (•faae)(x*-dx) ) ; 


l oca.  Varlab.aa 

,  cBxyl.x-fn.)  "get 

■ad 


-O  -o  eoaap.pa.Btt  -g  coaqj.pa.Bit . c  -I/aar/taclada/.ocal  -lp*atllb  -11*  -1*" 

-c  coop. pa. Bit  -g  co*p_pa.Blt . e  -I/aar/laclada/local  -L/opt/llb  -Ipaatllb  -It*  -1*“ 
-o  dco*p.pa.Blt  -g  co*p_pa_Blt .c  -I/mar/lacloda/ local  -L/opt/llb  -Ipaatllb  -11*  -1*“ 


Appendix  C 


Programs  Used  to  Evaluate  the  Matched  Filter  Bound 


Algorithm  2  Compute  data  for  matched  filter  bound 


•  CVat  lo 

•  :ta!  Hata  2at*ara1ty 

•  Data  Oct  jm .  i lew 


. . . / 

■lac.ada  <p*at . te  fe> 
lltt.tda  <*. lb  »> 

•lac.ada  <aaadtypaa  b> 

■lac.ada  (atdlo  h> 

•lac.ada  <«td.lb  b> 

•lac.ada  Caatk  »» 

I  lac.ada  Catrlmg  b> 

•  lac  Lada  tcoap.x  X> 

■laclada  <gTa*d  X> 

•lac.ada  <arati.. k> 

■daflaa  *la(a,b)  Ut)  <  (b)  *  Cx)  :  (b)) 

Bdaflxa  aax(a.b)  Cx)  <  Cb)  ’  (b)  :  (x)) 

•daflaa  ealloc.array  pcallot. array 
•daflaa  cal loc. Matrix  pea. ioc. Matrix 
•daflaa  DX.LO  26000000 
■daflaa  M.prob.atap  0.001 
■daflaa  r2  1  41421 J6B2J7 

/•  global  aarlab’.aa  •/ 

lat  typa;  /•  1:  alagla  fametloa,  2:  two  faactloaa.  •/ 
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anto.lnngth; 

croti.lnogthl ; 


lot  flnd.nv; 

donbla  *aoto_tnbln,  anto.nnlt,  auto.ttnp, 
lot  anto.no.tnrmf ; 

donbla  •ero««_tnblnl ,  crott.nnltl,  eroit_«tapl, 
int  cro«*_no_tnrMi 1 ; 

donbla  *crott_.tnbln2,  crc««.nalt2,  croaa.atap2,  croi*_lnngth2; 
lnt  tro*t_no_tarM«2; 

Int  no.pntb;  /*  nnmbar  of  path  •/ 

donbla  *c_tnn,  *alpba:  /»patb  dalay,  path  root-mann-tqi  magnltnda  •/ 

donbla  /•  h  th*  auto.corr  Matrix;  I  tha  croaa.corr  Matrix  •  / 

donbla  ■ymbol.int.rn;  /•  rnlativn  valua  (unit)  of  tha  tlgnal  •/ 

donbla  tymbol.lntarval ;  /•  abaolnta  ralnn  (unit)  of  tba  tlgnal  •/ 

dcoMplz  *Z; 

donbla  xbs,  zlx; 

donbla  *dkO,  *dkl; 

donbla  •dJrt.prob,  *<Htl_prob; 

donbla  DK.prob.rnnga; 

lnt  di.prob.lnngtb ; 

donbla  SJ; 

donbla  latagratlon.ranga; 

/ . / 

void  lnlt_antocorr_tabla(chnr  •); 
void  ialt.crottcorr. tabla (char  *,char  •); 
void  lnlt.channal (cbar  •)  ; 
donbla  fix  (donbla); 


donbla  Rry (donbla); 
void  Inlt.H(vold) ; 
void  lalt.l(vold) ; 
void  lalt.2( void) ; 
void  rand.Z(vold) ; 
donbla  cal. xMx( void ) ; 
donbla  cal. xlx( void) ; 
void  cal. d*0( void) ; 
void  ca..dbO_dlrl  ( void) ; 
void  hpaort (nntignad  long,  donbla 
doab la  Q ! doob la) ; 
donb.a  qf (donbla); 
void  tilt _0K _prob (void) ; 
doab.  a  P* .dbO ( doobl a!  ; 
doab. a  PZ.dkUdOnb.a)  ; 
doab. a  fnac .prob.dbO (donbla) ; 
doab. a  fnat. prob.dbl (donbla) ; 
told  ontpntlcbar  •.  cbar  •); 

. . / 


/•lnltlata  aoto.comlatloa  tabla*/ 
/•lnltlata  croat.cormlntlon  tabla*/ 
/•lnltlata  cbaanal  ckarn •/ 

/•anto.corr  nitk  mnlty  aalt  /sjM_lnt  •  / 
/•cront.corr  vltk  unity  nnlt  /  tym.lnt  •/ 


/•ganarataa  coMplaz  Gan  it  Ian  randon  »  In  Z  »/ 


doab. a  • lap. a (donb.a  ( *f nnc ) (donbla) »  donbla,  donbla); 
doab . n  da . tat ; 

. . / 


mata(1at  argc.  cbar  «*arg*H 

cbar  •antoeorrf  sin*  .  *croatcorrf  liaal ,  *cro«»corrf  nama2  .  *ehannalfaaaa,  *ontpntfnaM ; 

lat  iitteomuMClar ,  croiicormtatcktrl ,  crotacorrnaaacbar2,  channalaamacbar ,  oatpntf  aasacbar 

1* large  —Hi 

pr  1  at ' ( "aaaga  eoMp.pn  paravf  1  laaamaSn**)  ; 
atit 'o: . 


atraad  ar  gv [  i  J . “aym. lnt ” , k»y»6o 1 . latarval . DOUBLE  1 ; 

titseorr'iMH  •  (cbar  *JavTaad (argv [ 1] , “antocorr" .kantocorrnaaacbar .CHAR) ; 

crotacorrf  aamal  •  (cbar  • ) avraad  (argvtl ]  ."croticorr  1  “  .berot  acomiamncbnr  1  .CHAR)  ;  /•  •  nlda  •/ 

troatcorr*maa»a3  ■  (char  •)  avraad  (argv[  Ij  ,*’cro«icorr2'' ,hcroaacorrna»achar2.CHAJt)  ;  /•  -  (Ida  •/ 

ontpnt'ia—  ■  (cbar  • )  avraad  (argr[  1  ]  ,**oatpntf  Ha"  .boutpntf  namaebar ,  CHAR) ; 

cbaaaa . fiaaa  •  (cbar  •) avraad ( argr( i ] .  "cbaaaa 1 " .bcbannalnamacbar .CHAR) ; 

atraad  argt [ i J . "Da . t a*" ,bda . tat ,OOURUO ;  /•  lntagratloa  ttnp  •/ 

atraad  argt  [  1  J  .  “t  fpa"  .*t  ypa  .  (XT  )  ; 

atraad  argt { 1 J .“find. at" ,tf lad. at , 117) ; 

i a 1 1 .antocorr  . tab .a . avtocorrf tiat) ; 

Iftypa  —  3 

tyabo . . lat _rn  •  2*tyW>o . . latarva. ; 

a.  »a 

ty—o  . .  tat . ra  ■  aywfeo.  .latarval ; 
if(typa  -  a; 

talt . crottcorr  tab crotacorrf aa*a I . crotacorrf aa—2) ; 
la  It .  cbaaaa  .  '  cbaaaa .  *»aata  .  . 
lalt .n :  ; 

IfUypa  “  2. 

talt .1 C  1  . 
lalt. 1(1. 

Ifttypa  —  I  : 
c*..dto: ) . 

a.  aa 

ca..db0.dbt : ; . 
lalt.ox.prob ( i , 

oat  pat  ( oat  pat  fan—,  cbaaaa  1  fan—)  ; 

► 

void  lalt. antocorr. tab. a(cbar  *fliaaama) 

{ 

r:m  ♦tafiia, 
donb.a  taa . 
lat  1 , 

/«  antocorr  ft. a  format: 

ant: 

total.',  nagth 
•tap. alia 
no. of .tarmt 
tan  va.aa 


•/ 

Inflla  ■  fopaa(f liana— .“r") ; 
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if  (iaflla  «  IUUL)  < 

printf (“Cottld  sot  op*n  aatoeorrf ila\n") ; 
•lit (0) j 


f a caaf  (Iaflla, "Xlf \a” ,**nt o_aa it ) ;  /*  b«le  «ait  langtb  •/ 

dcuf  (lnfil*,"Xlf\n",fcanto_laagth) ;  /•  •  of  aait  tha  fil«  coataiat  •/ 

ficuf  Ciaf  lla/'XlfW.Eauto.ttap)  ;  /•  ttap  *i*a  la  aait  batvaaa  tarma 

f  a  canf  ( iaf  il« ,  “Xd\a  ” .  tanto.ao. tarma ) ;  /*  »  of  t«nu  •/ 

auto.tabla  -  (doobla  •  ) calloc.arrayC  DOUBLE, anto_no_t arma ."ant o.tabla") ; 

f  or  ( i«0 ;  Kanto.ao.tarm* ;  1+-0  { 

f acaaf  (iBfila,"XlfXlf\n",*tan,  Bauto_tabla[i] ) ; 


> 

f clota(laflla) ; 


> 


*/ 


void  talt.croticorr.tablaCchar  •filaBamal,  char  •filanama2) 

{ 

FILE  « inf 11a: 
doabla  tan; 
lnt  1; 


/•  crottcorr  fila  format; 
•alt 

tot *1.1 anyth 
lt«p.»lM 
ao_of .tarmt 
tan  vala* 


•/ 

laf  11a  ■  fopaalf  ilaaamal  .“r")  ; 
if  ( laf  l  la  “  rtILL)  { 

prtatf  (“Conld  aot  opaa  croaacorrl  fila\aM); 

•ilt(0) ; 

> 

f  acaaf  (laf  11a. “Xl f \a''  .tcrota.aalt 1) ; 
f acaaf ( laf  1  la  ,"Xif\a‘' .Ecroaa.laagthl ) ; 
f acaaf (Iaflla . “XlfNa" .bcroat.atapl ) ; 

♦acaaf ( iaflla, "Xd\a“. tcroaa.ao.tarmal) ; 

crota. tab. a!  ■  (doobla  •)calloc_array(!XroBLEpcrot«_BO_tarMl,"cro«a_tablal"> ; 
for  ( l«0 .  l<croat  .ao.taraal ;  1»*H 

f acaaf ( laf 11a. “XlfXIf \a” .fctao .  »croaa_tablal [ij ) ; 

! 

f c.ota! laf 1  la) ; 

lafl.a  »  fopaa  Cf  l.ao*a»a2  .“r”)  ; 

If  1  laf  1  .a  —  HU)  { 

prtatf (“Coaid  aot  opaa  eroaacorrt  f ila\i") : 

•alt (0) . 


f acaaf ( laf t .a.“Xlf\a” ,icro*a.aalt2) ; 

f acaaf (lafl.a. "Xlf \aM.*croaa_iaagth2) ; 
f acaaf (laf 1 la. “Xlf \a“ ,»croaa_atap2) ; 
f tcaaf (laf l.a.“Xd\a".ncroai_ao_tarma2) ; 

croaa . tab . *2  •  (doab.a  • ) cal  1 oc .array ( DOUBLE , croa a _ao  _tarm«2 , *' Croat _t ab 1 *2” ) ; 
for! t  *0 . I«croaa_ao_tarma2:l»*){ 

f acaaf ( laf l .a . “X. fXlf\a“ ,ttaa.  *croia_tabla2[lJ ) I 

fc.oaa! lafl .•; . 


void  lalt.cbaaaa. (char 


•cbaaaal ) 


F:l»  •lafl.a. 


/•  cbaaaa.  fl.a  format 
ao.patb 

c.taal  a.ptal 
C.taa2  a. pa* 2 


lafl.a  •  f opaa ! caaaaa.  , “r“ ) . 

If ( laf  1  la  —  r-XL)  { 

prtatf ( “Coa .d  aot  opaa  cbaaaa 1  fllala*); 

•tit (0) , 

> 

f acaaf (lafl.a  .“Xd\a“ .tmo.path) ;  /•  aombar  of  patkat  •/ 

c.taa  ■  (doab.a  • »c*. .oc. array (  XKJbl-I  , ao.patb ,"c_t*a") ; 
alpha  •  (doab.a  •! ca. .oc. array (DOUBLE . ao.patb , "alpha") ; 
for  (1M>:  t  <ao.pat b  ,  I**) 

f  acaaf  ( laf  1  la  ,“X.  fX„  f  \a“  .bc.taalt  J  .  (UIpbalO)  ; 
fcloaadaf  lla)  . 


told  la  It. Z (void) { 

Z  ■  (dcompli  *)mal loc (tlaaof (dcompli) “ao.patb) ; 


▼old  raa J.Z( void) ( 
lit  1; 

for ( 1 *0 ; l<ao_path ;1**5 

Z[l]  ■  dComplit  (grab  ()/r2)  .(graaO/rt)) ; 

> 


doabla  Rx (doabla  tat) 

i 

doabla  t.  tcalad.taa,  dalta; 
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lot  n; 


tan  /■  ayabol.lnt.ra; 

1# (t«u>0) 

t  ■  tan; 

alaa 

t  »  -tan; 
t  •«  anto_nnit; 
acalad.tan  ■  t  /  anto_*tap; 
n  »  (lnt)  acalad.tan: 
if  (n>*<nato_no_t*raa-0 ) 

r*tum(anto_tabl*[anto_no_tara»-l] )  ; 
ala*  { 

dalta  ■  acalad.tan  -n; 

ratnrn (anto.tabl* [nj  ♦  dalta* (anto_tabla[n-*t]-auto.tabla[a] )) ; 

} 

> 


donbla  Rxy (donbla  tan)  /*Croaacorr*latlon*/ 

{ 

donbla  »croaa_t*bla ,  croaa.ult,  croaa.atap; 

tnt  croai.no.taraa ; 

donbla  t,  acalad.tan,  dalta; 

lnt  a; 

tan  /■  ajabol _ lnt.ra ; 

lf{tnnX>H 

t  ■  tan ; 

croaa.tabla  *  croaa.t^blal ; 
croaa.nnlt  •  croaa. nnltl; 
croaa.atap  ■  croaa. atapl; 

Croat. no. taraa  a  croaa.no. tana  1 ; 

> 

*:**{ 

t  «  “tan; 

croaa.tabla  -  croaa. tabla2; 
croaa.nnlt  •  croaa. nnlt2; 
croaa.atap  •  croaa. tt*p2; 

Croat. ao. taraa  ■  Croat _ao_taraa2 ; 

> 

t  ■ *  croaa.aalt  , 

tea. ad. tan  at/  croaa.atap; 

a  •  (tat)  acalad.tan; 

t f  J n>- ( croaa . no . taraa - 1 ) ) 

rat  nr* (croaa. tab ;*[croi a. no. taraa- l] ) ; 
alaa  { 

data  ■  acalad.tan  -a; 

rat  nr*  (croaa.  tab  l*t*J  *  d*lt**(croaa.tabl*[n-*l]-croaa_tabl*[a])) 


void  lalt.K (void? 

( 

tat  l.  J. 

doab.a  ■  •Qa,  **L*abda,  •*.*,  ana; 

n  *  (doab.a  « • ) cal loc.aatrlx C30CBLK , no. path , no. path , "H") ; 
for (t»0, 1< no. path; 1**J 
for! j-0 ;  Jtao.patt ;  )«-») 

"UJ  [  jj«*.pa*llj*alpaaljj  *Rx(  (c_ta*[lj-c.tantj])) ; 
print? (“a . \a“!  . 
f  or  !  1  "0  ;  1  <*o.p*ta  ;  !*•*){ 
f  or!  J*>0.  J  cao.pat h  ;  J**  ) 
prtatf CXa  ".«lljtjj> . 
prtatf !  , 


•old  lalt.l (rotdl 

tat  I,  j. 

1  •  (doab.a  ••ica..oc.Mtrli(  0O1TBI.J  ,ao. path. no. pat*  ,  ”¥")  ; 
for  (1*0,  t<ao_ pat*  .1*-*; 
for  ()»0;  J<ao  .pat*  ,  J 

*l»J  l  Ji"«  .p*aU)  *a.p*a(  J)*Rt7(  ( c. tan lij -c. tarn [JJ ) ) ; 

> 

doab.a  ca. .a*a ' ) 

( 

tat  1.  j. 

dcoap.t  taap,  taapl; 

taap  1  ■  dCoop.xtO,  0), 

for (1*0;  ICao.pach;  1*-*) 
for ( J*0 ; jtao.patb ;  J«~*)  i 

taap  -  dCooplxt *lij  t  J  J  .  0); 
t ootp  ■  desalt taap.  Z[tJ); 
taap  •  dConKtaap.  dCoaJg  ( Z  l  JJ  ) ) ; 
taaipt  ■  dCadd  '  taap  ,  t«atpl), 

> 

ratnr* (taopi  ra) ; 

> 

donbla  cal.alaO 

< 

lnt  1.  J; 

dcoatplx  taap,  taatpl; 


taap  1  ■  dCoapl*(0,  0)  ; 


for  ( 1*0 ;  1< no. path  ;!«-*} 


f or ( J -0 ; j <»o_p*t h : J ♦♦ H 

taatp  •  dCoaplx(J[l]  [j]  ,  0); 
taap  ■  dCkraKtoaqj,  Z[i] ) ; 
tap  •  dCanlCtaap,  dCon jg(Z[j] ) ) ; 
tap  1  ■  dC*dd(t*ap,  taapl); 

> 

rttnrt (taapl .ra) ; 

> 

▼old  c*l_dk0_dkl  O 

i 

lot  ft; 

donbl*  aqr.xK*,  t,  b  ; 

dkO  ■  (doabla  •)c*lloc_*rr»j([XnJ8LB,  (DK_LHT*1)  ,“dlr0“); 
dkl  ■  (doabla  •  )t*lloe_*n,*j(DOUBLS,  (DK.LD+1)  ,*'dklM) ; 

f or ( 1 -0 ; 1 <»DK _ LEI ; i+*) { 
rud_Z() ; 
x*m  »  c*i_*kx()  : 
ill  -  c*I_*l*C); 
aqr_*Jl*  ■  aqrt(xJtx); 
dkO[l]  ■  ( x.Hx  -  xls)/aqr_xflx; 
dkl[t]  ■  (*Ax  •*  **x)/»qr_x«x; 

/•  Ihsm  1  It  aont  «t  thla  ica1«  ahiia  kO:  -1  on  otfcar  ac*lo 
kl :  1  ob  othar  ical* 

•/ 


> 

hptortOK.LO,  dkO)  ; 

kpiort !  DK  _LU ,  dkl): 

for(l-0;l<DK_LKI;l~M 
dkOtlJ  -<lkO[l*lJ  ; 
dklltj»dkl[l«l]; 

> 

•  -  dk0[  (Dk.LCT-1 )  ]  : 
t  •  dkl  [(Dk.LU-l)J  ; 
tatograt to*. rang*  •  Hi(i,b); 

3*  .  prob .  ru(«  ■  1  atagrat  Ion.  rug* ; 


»«td  c*..dkO;j 


d*0  •  (doab.o  •  Jcalioc.arrey  (DOUBLE,  (DK.UaM)  ,HdkOH) ; 
for:  t-0. l<»Dk_ ld :!•*-»)( 
raod.ZC  . 

•w*  ■  t* ; .  l«i  ( ) ; 

dkC[  t  J  •  iqn  («n*J : 

*paort  :Dk. LO.  dkO); 
fortl-O.  tOl.LD,  1~) 
dkOll  j-dk«lt*lj  ; 

latograt  too.  ruga  ■  dkOl(DK.LD-l)]  ; 

3k.pro0.rugo  ■  latograt  loa. ruga; 


▼old  talt . 3« .prob w 

tat  «M  prob.l.  dkl.prob.t,  dkO.t.dkl.l,  bogla_i,  oad.l; 


dk  .pro*,  .oogt a  •  '  tat ; 3k.prob.rugo/DK .prob.atop; 

11  profc  .  oogta  ao  |  . 

d*C.pro*  •  :»oob  o  ■  :  ca . .  ot  _ array : 3GUBUI , dk _ prob _  ioagt  h ,  "dkO.prob" ) ; 

do. to  ■  3«  pro*. at op. 

dkO.  t  »  0 
dkC_prot[0)  •  0. 

for  : itt  pro*,  t  ■  l  .  d*}.pro*.1<adk.prob_  ugtk  ;  dkO.prob.  1** ,d*lt*+BDK_prob_atopH 
bog  la  1  *  dkO. I . 

in  o:  :d*fl(dk0.tJ<do.t*.«'.(ik0.1<CDk.Ln-l)J ) 
dkO.t*-. 
ood  .  i  •  dkO  .  t , 

dkO.  pro*  1 1W.  pr  o*  .  1  j  •  ;  doo*  .  •  !  !  ood  .  t  -  b*gi  a_  1 )/ (DK.LO*  DK.prob.atop)  ; 

/ 

pfroo.arroy  d*«: . 
tftt^po  —  2) * 

dk  1  .pro*  ■  :  doob  .  a  •  ;  co.  .  oc  .array  !  3CCBL1  .dk.prob.  lagtl , “dk  t  .prob" )  ; 

do.ta  ■  3k .pro*. atop , 

dkt .1  bo. 

dk I .pro* [0 J  o  0. 

for (dk  l  .pro*,  t-l . dk  1  .pro* .  t  <adk  .prob.  1  ugtk :  dk  1  .prob.  1*-* ,  dal  ta*-DK_prob.  atop )  { 
bogla. t  •  dk 1 .t ; 

a*l  .o'  !dkl  (dkl  .tJ<4o.to)U!dkl.l<(DK_UB-l))) 
dkl  . 

OOd_ 1  *  dkl.l. 

dkl .probldkl.prob.l]  •  (doobl*) (**d_l  -  bogia_l)/(DK..UtI»Dk^prob_itop) ; 

> 

pfro*. array ! dk  t ) . 

> 


doabla  PF.dkOCdoablo  f  ) 

i 

lat  •; 

doablo  acalod.x; 

If  £r—0) 
rotara(O)  ; 

atal*d_*  ■  t /Dk.prob. atop ; 
a  ■  (lat)  ac*l*d_»; 


If  (»>B(dk_prob.lugt*-l) ) 


return (0) ; 

•Is* 

retura(dkO_prob[n-*l]} ; 


doable  PB.dk 1 (double  z) 

{ 

lit  ■; 

double  acaled.x; 

If  (i»0) 
return (0) ; 

scaled.*  »  x/DK.prob.step ; 
i  ■  Clat)  scaled.x; 

If  (n>«(dk_prob_leugth“l) ) 
return (0) ; 

else 

return  (dkt_prob[n+l]} ; 


double  fuac.prob.dkO (double  E) 

{ 

double  10.  teap ; 

leap  -  -(SI/10) ; 

10  ■  pou( lO.taap) ; 

teaip  •  sqrtCIO);  /»Eb«lt  B**»2  «/ 

teap  ■  Q((E/teap))«PB_dkO(B) ; 
return (leap) ; 


doub.e  fuac.prob.dkl (doable  E ) 

{ 

doub.e  V0.  teap; 

teep  •  ‘(SI/101 ; 

■0  •  poe! 10. teap) ; 

teap  •  uqrt(IO);  /»Eb-l,  Is<2  */ 

teap  -  0 ( ! K/teap) )"Pt_dkt(S); 
return  teap  )  ; 


sold  output (cber  *fneae,  ebar  •cfneae) 
KLI  »fpt. 

tier  output  ft .e (30) ; 
doub.u  Integra. ,  ever age IbOO] ; 
let  I  . 
doub.e  tee. 


•  prtstf : out put f I . e.“Xs . Xd-X . teXs  pen— , f aeae.tfpe , symbol. Interval , cfui 
fpt  •  f  opee  :  out  put  f  1.  a  ,***“) ; 
if (fpt  —  r-lL  1  « 

pr t at f ' “Con . d  eot  opea  output  file!  \n"}; 
ettt : o; . 


ik>. 

for  SS-c.  SB<2i>,SS— 0  .  1)( 

Integra.  ■  *  lap .e ( fuse .prob.dkO ,0. integretlon.range) ; 
If.trpe  —  i ; < 

fprtatf : fpt . ”Xf  Xe\e" .SI, lstegrel ) ; 
ff  .u»S  ’ ivu  ,  . 

e*erege(i*«J  •  letegre.. 


If :  tfpe—  2:  • 

1  ■  0. 

for  ;S»-b  ,  M<2i..  SI-eO  t ) «. 

Integra  •  s tap . e(f eec .prob.dkl ,0, let egret  lorn. range} ; 
e*eragaH~*J  ••  Integra. . 

I  •  0, 

for  :  SJ-t» .  ,  ss**o  lit 

fprtstf ; fpt .  ”Tf  1*\*~ .SI , ; average ( l**J/5)!; 


f e . ose ( fpt ; . 

uprtetf  (output fl .»,“!«  Id-I  lels  .  pdf“  .fnnae.t n>«.»r»»>°i..latarTel  .eft 
rpt  •  fopeeioetputf t ;e,“u“) ; 
if  (fpt  —  mi:  * 

prlutf (“Cou.d  eot  opes  output  file!  \»*}; 
etlt (0) . 

> 

fprtstf (fpt , "Tit leTert :  POP  Spa. Internal :Xf\s" .symbol. Interval) ; 
Integra:  •  a  lap: e (Pl.dkO .0 „ ietegrst lou.raage} ; 
fprlntf  Cfpt.“Vdk0  1  Xf\"  \e", lstegrel ); 
for (taooO;tas<etategrattoe.raage;tau  *e  0.1} 
fprletf (fpt .“If  If \a" .tee .Pl.dkO (tan) } ; 

If  (typo— 3)  < 

Integral  ■  ttap. elPI.dk t ,0, lu t egret tos.raage) ; 
fprtstf (fpt ."\a\“dk t  *:Xf\“  \a" , Integral ) ; 

f or (teu-0 . teu<«t»tegruttoo_reage; tea  0.1} 
fprtutf (fpt .“Xf  If \n"  .ten  .PI.dk 1 (tas)) ; 

> 

fclose(fpt) : 

> 

/ . 


Numerical  Integration 


r 


. . . . . . »• . . . 

donbla  alapla (doobla  (*fnnc) (donbla) ,  donbla  a,  donbla  b) 

{ 

doabla  tao,iotngral»0; 

for  ( t  mm  ;  t  an<»b ;  t  na-»«dal  t  at ) 
intagral  *-  (‘fane) (tan) ; 

integral  daltat; 
ratnns  intagral; 


donbla  CKdoabla  *) 

i 

donbla  a.b.e; 
b*r***£ .61 ; 

a»0.86«iqrt(b)*l . 66»i; 
c« (donbla) (-i*i/2) ; 
t» (donbla) axp(c) /a; 
ratnm(i)  ; 


void  fcpaort (analgnad  long  a,  donbla  raU) 

i 

aaalgaad  long  l.lr.j.l; 
donbia  rra; 

If  (a  <  2)  ratara; 

>!a  »  l)*l ; 
ir»n : 

for  ' ; ; )  { 

If  l.  >  1)  { 

rra»ra[--l] ; 

*  a.aa  i 
rr»»r*Llr J ; 
raltrj-rallj . 
if  tr  —  I)  i 
ral i] *rra 


aftt.a  '.)  <■  ID  * 

If  ]  Mf  U  ratjj  <  ra l J*  1  ]  5  J-*-*; 
if  '.rra  <  raljjl  * 
ra(ij*ral )} . 

j  <«-  l, 

•  a.aa  )«tr-l . 

rail J  »rr« 


iwaa-trj  t«ariu.  ftaclpaa  Software  m.l.  •/ 


Lof.a .  Vartan  .  a« 
co^l.*-'.«M  aptc  -O 
ena^t  •■».(«•»«  "apct  -0 
'gtt  -0 


-g  toap_pa . e  -I/aar/taclodo/local  -I/hoan/ilJrap/boanl/intlodn  -L/opt/llb  -Ipaatllb  -11a 
-g  comp  pa . c  - I /asr/laclad a/local  -L/opt/llb  -Ipaatllb  -1 la  - la  -lcoapir'* 

-g  toap.pa.t  -I/nir/lnclndn/locol  -I/boWainap/boaet/lnclnde  -L/opt/llb  -Ipaatllb  -11a 


-1*  -Icoatplr" 
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